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Abstract: We analyze all available solar and related reactor neutrino experiments, 
as well as simulated future 7 Be, p — p, and pep solar neutrino experiments. We treat 
all solar neutrino fluxes as free parameters subject to the condition that the total 
luminosity represented by the neutrinos equals the observed solar luminosity (the 
'luminosity constraint'). Existing experiments plus the luminosity constraint show 
that the p — p solar neutrino flux is 1.02 ± 0.02 (la) times the flux predicted by the 
BP00 standard solar model; the 7 Be neutrino flux is 0.93+063 the predicted flux; 
and the 8 B flux is 1.01 ± 0.04 the predicted flux. The CNO fluxes are very poorly 
determined. The neutrino oscillation parameters are: Am 2 = 7.3lg;g x 10~ 5 eV 2 and 
tan 2 6*i2 = 0.41 ± 0.04. We evaluate how accurate future experiments must be to 
determine more precisely neutrino oscillation parameters and solar neutrino fluxes, 
and to elucidate the transition from vacuum-dominated to matter-dominated oscil- 
lations at low energies. A future 7 Be v — e scattering experiment accurate to ±10% 
can reduce the uncertainty in the experimentally determined 7 Be neutrino flux by 
a factor of four and the uncertainty in the p — p neutrino flux by a factor of 2.5 
(to ±0.8%). A future p — p experiment must be accurate to better than ±3% to 
shrink the uncertainty in tan 2 Qyi by more than 15%. The idea that the Sun shines 
because of nuclear fusion reactions can be tested accurately by comparing the ob- 
served photon luminosity of the Sun with the luminosity inferred from measurements 
of solar neutrino fluxes. Based upon quantitative analyses of present and simulated 
future experiments, we answer the question: Why perform low-energy solar neutrino 
experiments? 
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1. Introduction 

We assess in this paper how well existing and future experiments can determine 
neutrino parameters and solar neutrino fluxes. The starting point for our calculations 
is the assumption that the large mixing angle (LMA) MSW solution |T| describes 
exactly the behavior of solar neutrino oscillations and that the true values for the 
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oscillation parameters lie close to the current best-estimate oscillation parameters 0, 

a 1 - 

The conclusions derived in this paper are based upon 152 global analyses of all 
of the available solar neutrino and reactor data. The individual analyses differ, for 
example, in the input data (including cases before and after the recent announcement 
of the SNO salt-phase data and the improved GNO and SAGE measurements), the 
different constraints placed on the neutrino fluxes (e.g., with or without the luminos- 
ity constraint), the number of neutrino fluxes that are treated as free variables, and 
the different assumptions made about future solar neutrino and reactor experiments. 
The purpose of this introduction is to set the stage for the many detailed results that 
are presented in the subsequent sections and to provide the reader with a overall map 
of what is contained where in this paper. 

The new experiments discussed in this paper will test whether the LMA neutrino 
oscillation solution is correct and, if it is correct, will determine accurately the solar 
neutrino fluxes, the solar neutrino mixing angle, and the difference of the squares of 
the neutrino masses. We will not be surprised if these experiments reveal physics or 
astronomy that cannot be explained within the now conventional framework of the 
LMA oscillation solution and the standard solar model. The fundamental goal of our 
paper is to determine what is expected in order to make it easier to recognize what 
is unexpected. 

We also show that neutrino experiments alone can be used to measure the current 
rate of nuclear energy generation in the Sun and to compare this neutrino-measured 
solar luminosity with the conventional photon-measured solar luminosity. This com- 
parison will provide an accurate test of the fundamental assumption that nuclear 
fusion among light elements is responsible for the solar luminosity. Moreover, this 
same comparison will test a basic result of the standard solar model, namely, that the 
Sun is in a quasi-steady state in which the current energy generation in the interior 
equals the current luminosity at the solar surface. The heat being radiated today 
from the solar surface was created in the solar interior about 40,000 years ago j|. 

Although most of our effort in this paper is directed toward new experiments, 
the most surprising result is that the p — p solar neutrino flux is determined to ±2% 

lr The data included in the analysis in the present paper include all relevant solar and reactor 
neutrino data available on or before September 7, 2003, the date on which the SNO collaboration 
announced the results of their measurements in the salt phase and on which the GNO and SAGE 
collaborations announced refinements of their measurements. The recent SNO, GNO, and GALLEX 
data, which are included in the analyses reported in this paper, were all presented at the VHIth 
International Conference on Topics in Astroparticle and Underground Physics (TAUP03), Seattle 
(Sept. 5-9, 2003). Only minor numerical changes resulted from this updating, but we are grateful 
to the referee and to JHEP for allowing us to include the new data in the final version of this paper. 
An energetic reader can learn interesting and encouraging details about the robustness of current 
solar neutrino inferences by comparing the many numerical results presented in hep-ph/0305159,v2. 
with the numerical results presented in the present version of this paper. 
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by the existing solar and reactor experiments if one imposes the constraint that the 
luminosity of the Sun is produced by fusion reactions among light elements that also 
produce solar neutrinos (the so-called 'luminosity constraint'). The best-fit value of 
the p — p flux inferred from a global fit to existing neutrino experiments and the 
luminosity constraint is within 2% of the flux predicted by the BPOO standard solar 
model. 

1.1 Historical perspective 

The first forty years of solar neutrino research has demonstrated that new physics 
may appear when we carry out neutrino experiments in a new domain of sensitivity. 
Most of the new experiments considered in this paper will be sensitive to neutrino 
energies that are less than or of the order of 1 MeV, a domain in which solar neutrino 
energies could not previously be measured. More than 98% of the predicted flux of 
solar neutrinos lies below 1 MeV. 

One of the primary reasons for carrying out the new experiments is to test 
whether new physics shows up at lower energies. As stated earlier, our explorations 
in this paper are intended as a guide to the expected in order that they can be used 
to help identify the unexpected. 

Moreover, we want to use solar neutrino experiments for their original pur- 
pose 0: to test theoretical models of how a main sequence star gains energy and 
evolves by burning hydrogen to helium. To do this, we must have sufficient experi- 
mental data to measure the total neutrino fluxes that are produced by the principal 
neutrino-producing nuclear reactions. Until very recently, it was necessary to assume 
the standard solar model predictions for all the solar neutrino fluxes and their uncer- 
tainties in order to determine reasonably constrained values for neutrino oscillation 
parameters. Only after the results of the Super-Kamiokande and SNO experiments 
became available was it possible to measure directly the 8 B neutrino flux. However, 
even now, the other solar neutrino fluxes must be taken from standard solar model 
calculations in order to obtain the well-constrained global solutions for neutrino os- 
cillation parameters that are seen so frequently in the literature 0, |3|. 

A study similar to what we carry out in this paper was performed in 1996; the 
earlier investigation provides historical perspective. In 1996, only the chlorine [0], 
SAGE [0], GALLEX || and Kamiokande |J experimental results were available. In 
a paper entitled [ID 'How well do we (and will we) know solar neutrino fluxes and 



oscillation parameters?', the authors found seven years ago that the small mixing 
angle (SMA) solution was slightly preferred over the LMA and vacuum oscillation 
solutions. All three, SMA, LMA, and vacuum oscillations, were allowed at la. The 
uncertainties in the solar neutrino fluxes were large. At 95% C.L., the uncertainties 
were a factor of two for the p — p neutrino flux and a factor of five for the 8 B neutrino 
flux. The 7 Be solar neutrino flux could be as large as 6.4 times the 1995 standard solar 
model prediction. In this investigation, there were only enough experiments to permit 
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treating one flux at a time as a free parameter. The other neutrino fluxes and their 
uncertainties were required to be consistent with the 1995 standard solar model [|ll| . 
The principal constraints on the p — p and 7 Be neutrino fluxes were established 
by the luminosity constraint [12], not the then existing solar neutrino experiments. 



In 2001, Garzelli and Giunti JT3J performed a Bayesian analysis independent from 
the Standard Solar Model (SSM) predictions for the solar fluxes. The luminosity 
constraint was used to obtain the posterior distributions oip—p, 7 Be and 8 B neutrino 
fluxes in agreement with BP00 predictions [[14]]. At 90% C.L., the allowed ranges were 



[0.99,1.08], [0.02,1.15], [0.62,1.22] for p - p, 7 Be, and 8 B neutrino fluxes (in units of 
BP00), respectively. 



With the extensive results of three additional experiments, Super-Kamiokande [15 



SNO [ITTJ, [17]], and KamLAND flq| , we now know-or think we know-a great deal more. 
But, the enormous change in the consensus view about solar neutrinos that has been 
caused by the data from just three new experiments provides an additional caution 
that surprises may appear in the future. 

1.2 The theme of this paper 

The theme of this paper can be stated simply: the astrophysics of the solar interior 
and the physics of neutrino propagation both deserve to be explored independently, 
not assumed. We investigate the extent to which existing data and future experi- 
ments constrain separately both the stellar astronomy and the neutrino physics. 

A key step in our analysis is the imposition of the luminosity constraint |I2"| , |T9| ], 
which implements in a global way for the Sun the constraint of conservation of 
energy for nuclear fusion among light elements. Each neutrino flux is associated 
with a specific amount of energy released to the star and therefore a particular 
linear combination of the solar neutrino fluxes is equal to the solar luminosity (in 
appropriate units). One can write the luminosity constraint as 

= <L1) 

where L is the solar luminosity measured at the earth's surface, 1 A.U. is the average 
earth-sun distance, and the coefficient ojj is the amount of energy provided to the 
star by nuclear fusion reactions associated with each of the important solar neutrino 
fluxes, The coefficients are calculated accurately in ref. |T2 . 



In this paper, we present results that are obtained with and without imposing 
the luminosity constraint, in order to illustrate the power of the constraint. 

We include in our global analyses reactor [|18[] and solar neutrino data f6|, |7|, jB[ [20 



15, 16 , 17] . We marginalize all the results quoted in this paper over the mixing angle 
#13 using the calculated dependence of the global x 2 as a function of the dominant 
neutrino oscillation parameters, Am 2 and #12, the solar neutrino fluxes, and #13. For 
earlier calculations of the dependence of x 2 on #i3 ; see re fs- [0; 0) 0; especially 
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figure 7 of ref. [22| and figure 4 of ref. [23]. We use in this paper the dependence 



shown in figure || that is adapted from ref. [24 



We investigate in this paper the extent to which the inferred values of the neu- 
trino oscillation parameters are dependent upon the common assumption that the 
solar neutrino fluxes have the best-estimate values and uncertainties determined from 



the standard solar model [[TJj]. Throughout this paper, we present results in which 
the 8 B solar neutrino flux is treated as a free parameter. In the latter sections of 
the paper, we present results in which also the 7 Be solar neutrino flux, and later the 
p — p and CNO neutrino fluxes, are treated as free parameters. 

We begin our investigation by determining how well the existing set of solar and 
reactor neutrino experiments constrain the solar neutrino fluxes and the solar neu- 
trino oscillation parameters. Most of this paper is, however, devoted to experiments 
that are currently being planned or developed |^3J . We explore how well future solar 
neutrino experiments can be expected to determine neutrino fluxes and oscillation 
parameters, all provided there are no major new surprises. 



1.3 Organization of this paper 

1.3.1 Advice for generalists and specialists 

A quick glance at the table of contents will be useful for both generalists and spe- 
cialists. We have written the titles of the sections and the subsections so that one 
can easily see what subjects are covered and in what context. 

General readers will find everything they want to learn from this paper in our 
discussion section, section |9|. We answer in section |9.1| the question: Why do low 
energy solar neutrino experiments? And in section |9.2| , we summarize the principal 
results and ideas that are discussed in the main text. We give our view of what does 
it all mean in section |9.3| . 

There are many detailed results and explanations that could not be included in 
section |9|. For the reader who is a specialist in a particular area related to solar neu- 
trino research, we nevertheless recommend reading first the discussion section |] and 
only afterwards selecting a particular topic from the list of section topics described 
below or from the table of contents. 



1.3.2 Topics discussed in each section 

We present in section Q the approximate analytic behavior of the electron-neutrino 
survival probability, P ee , for the LMA oscillation solution. A number of the results 
in this paper can be understood qualitatively by referring to figure [I] of section |2|. 

Section | describes some of the technical details that we use later in the paper, 
including the input experimental data, the ingredients of the global x 2 > the rela- 
tionship between the two-neutrino and three neutrino survival probabilities, and the 
simulations we make of future KamLAND and SNO data. The details of section |3] 
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are primarily of interest to experts on neutrino oscillations and can be skipped by the 
general reader who just wants to know the scientific implications of our investigations. 

We use in section ^ all of the available solar and reactor neutrino data to deter- 
mine the current constraints on neutrino oscillation parameters and solar neutrino 
fluxes (see especially table f| and figure |3|). The global analyses in this section evolve 
in steps from treating only the 8 B neutrino flux as free, to letting both the 7 Be 
and 8 B fluxes be free parameters, to the most powerful method: treating all of the 
neutrino fluxes as free parameters subject to the luminosity constraint. We also an- 
alyze in this section simulated data for three years of operation of KamLAND and 
show how these additional data are expected to improve our knowledge of oscillations 
parameters and fluxes. 

We explain in section [| the reason why the luminosity constraint helps provide a 
powerful bound on the p — p solar neutrino flux, but imposes only a relatively weak 
constraint on the 7 Be solar neutrino flux. 

Section |5| and section [7| provide our best answers to the questions that exper- 
imentalists planning future solar neutrino measurements most often ask us. These 
sections address the parallel questions: What will we learn from a 7 Be solar neutrino 
experiment? What will we learn from a p — p solar neutrino experiment? In order to 
answer these questions, we analyze simulated data for 7 Be and p — p solar neutrino 
experiments. 

We show in section |6] that the current predictions for a 7 Be solar neutrino ex- 
periment depend sensitively on whether or not one assumes the correctness of the 
best-estimate neutrino fluxes, and their uncertainties, that are predicted by the stan- 
dard solar model. Table |B| shows that we will be stuck with a large experimental 
uncertainty in the flux of 7 Be solar neutrinos unless a 7 Be solar neutrino experiment 
accurate to ~ ±5% is performed. Fortunately, a 7 Be solar neutrino experiment will 
also provide a powerful constrain on the p — p solar solar flux via the luminosity 
constraint. 

We show in section [7] that a measurement of the p—p neutrino-electron scattering 
rate must be very accurate, ~ ±1%, in order to significantly improve our knowledge 
of the tan 2 9 or the solar neutrino fluxes (see especially table U). If this accuracy is 
achieved, then the allowed range for tan 2 8 12 will be reduced by a factor of two. We 
demonstrate also in this section that a measurement of the pep neutrino flux would 
be about as informative as a measurement of the p — p neutrino flux. 

We answer in section § two important questions: At what level is maximal mixing 
excluded by current data? How will the SNO salt-phase data affect the exclusion of 
maximal mixing? 

2. MSW at high energies, Vacuum at low energies 

The effective Hamiltonian for two-neutrino propagation in matter can be written 
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conveniently in the familiar form 



n, er m nn 



Am 2 cos26 1 2 



H 



V2G F n c 

4E 2 
Am 2 sin26i2 
2E 



Am 2 sin28±2 
2E 

Am 2 cos29i2 _i_ \/2Gf"c 
4E 2 



'2.1) 



Here Am 2 and 6 I 12 are, respectively, the difference in the squares of the masses of 
the two neutrinos and the vacuum mixing angle, E is the energy of the neutrino, 
Gp is the Fermi coupling constant, and n e is the electron number density at the 
position at which the propagating neutrino was produced. The best-fit values for 
Am 2 and tan 2 6 12 obtained from a global solution to all currently available solar 
neutrino and reactor anti-neutrino experimental data are (see table [|, the third row) 
Am 2 = 7.3 x HT 5 eV 2 and tan 2 9 12 = 0.41. 

The relative importance of the MSW matter term and the kinematic vacuum 
oscillation term in the Hamiltonian can be parameterized by the quantity, (3, which 



represents the ratio of matter to vacuum effects. From equation |2.1| we see that the 
appropriate ratio is 

= 2v^A 
Am 1 

The quantity (3 is the ratio between the oscillation length in matter and the oscillation 
length in vacuum. In convenient units, (3 can be written as 



(3 = 0.22 



1 MeV 



100 



g cm 



7 x 10- 5 eV 2 
Am 2 



[2.3) 



where p e is the electron mean molecular weight (p e ~ 0.5(1 +X), where X is the mass 
fraction of hydrogen) and p is the total density, both evaluated at the location where 
the neutrino is produced. For the electron density at the center of the standard solar 
model, (3 = 0.22 for E = lMeV and Am 2 = 7 x 10- 5 eV 2 . 

For the LMA region, the daytime survival probability can be written to a good 
approximation in the following simple form [27], |HT[ 



- + - cos 26^ cos 26 12 



where the mixing angle in matter is 



cos2# M 



cos 20i2 - (3 



12 



cos2# 12 -/3) 2 + sin 2 29 12 



[2A) 



(2.5) 



Figure |1| illustrates the energy dependence of the LMA survival probability, P ee . 
If (3 < cos 2#i2 ~ 0.4 (for solar neutrino oscillations), the survival probability corre- 
sponds to vacuum averaged oscillations, 



- sin 2 2#i2 {(3 < cos 26*12, vacuum). 



;2.6) 
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Figure 1: Schematic survival probability. The figure shows the electron neutrino 
survival probability, P e e; as a function of neutrino energy for the (daytime) LMA oscillation 
solution. For small values of the parameter f3 defined in equation |2.2| and equation ^7J, the 
kinematic (vacuum) oscillation effects are dominant. For values of (5 greater than unity 
the MSW (matter) oscillations are most important. For solar conditions, the transition 
between vacuum and matter oscillations occurs somewhere in the region of 2 MeV. 

If P > 1, the survival probability corresponds to matter dominated oscillations, 

P ee = sin 2 fl 12 09 >1, MSW). (2.7) 

The survival probability is approximately constant in either of the two limiting 
regimes, (3 < cos2#i2 and (3 > 1. The LMA solution exhibits strong energy de- 
pendence only in the transition region between the limiting regimes. 

At what neutrino energy does the transition take place between vacuum oscil- 
lations and matter oscillations? The answer to this question depends upon which 
neutrino source one discusses, since the fraction of the neutrino flux that is produced 
at a given radius (i.e., density and /x e ) differs from one neutrino source to another. 
The 8 B neutrinos are produced at much smaller radii (higher densities) than the p— p 
neutrinos; the 7 Be production profile is intermediate between the 8 Be and p — p neu- 
trinos. According to the BPOO solar model, the critical energy at which (3 = cos 2# 12 
is, for tan 2 12 = 0.41, 

E(crit) ~ 1.8 MeV ( 8 B); ~ 2.2 MeV ( 7 Be); ~ 3.3 MeV (p -p). (2.8) 
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The actual energies for p — p and 7 Be neutrinos are below the critical energy 
where they are produced. To a very good approximation, 8 B neutrinos are always 
in the MSW regime, while p — p and 7 Be neutrinos are in the vacuum regime. The 
energies are so low for p — p and Be neutrinos, that the energy dependence of 
the vacuum oscillations, oc sin 2 (Am 2 L / AE) , averages very close to one-half in any 
measurable energy interval. To a good approximation, the survival probability is 
given by equation pTO for both 7 Be and p — p neutrinos. 



The fact that the survival probability in equation |2.6| is independent of Am 2 
explains why the 7 Be and p — p experiments discussed in section ^ and section [7| are 
insensitive to Am 2 . We shall show in section ^ that a precise p — p experiment can 
be sensitive to tan 2 #i 2 . 

The rather low critical energy for 8 B neutrinos explains why spectral distortions 
are difficult to measure in Super-Kamiokande and SNO with present thresholds E e ~ 
5 MeV. However, a distortion should be measurable if the threshold can be lowered 
by 1-2 MeV. 

3. Some technical remarks 

In this section, we describe some technical aspects of our analysis and introduce 
some of the notation. This section can be skipped by readers who do not suffer from 
acute addiction to the details of neutrino oscillation analyses. Where necessary in 
the remainder of the text, we refer back to this section for definitions or for specific 
procedures. 

We summarize in section [371] (see especially table [l]) the existing experimental 
data used in latter sections to carry out global analyses of solar and reactor neutrino 
experiments. In section |3.2| , we describe the global x 2 we have used and define 
the free parameters in Xgiobai- m particular, we define in section ^72| the reduced 
neutrino fluxes, /b ( 8 B neutrinos), /ee ( 7 Be neutrinos), and / p _ p (p — p neutrinos). 
We define in section |3.3| the effective survival probability after marginalizing over 
# 13 . In section R75], we describe how we simulate KamLAND data for three years of 



operation. Finally, in section 3.6 we describe how we simulate the SNO data that 



will result from the salt-phase measurements. 
3.1 Experimental data 

Table [l] summarizes the solar (S) || [5], || [16, [17], |2(| and KamLAND (K) data 



used in our global analyses that are presented latter in this paper. The total number 
of data used in the analyses is 96; the number of data derived from each experiment 
are listed (in parentheses) in the second column of the table. In the third column, 
labeled Measured/SM, we list for each experiment the quantity Measured/SM, the 
measured total rate divided by the rate that is expected based upon the standard 
solar model and the standard model of electroweak interactions. 
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Table 1: Experimental data. We summarize the solar (S) and the KamLAND (K) 
data |l8}| used in our global analyses. Only experimental errors are included in the 
column labeled Result/SM. Here the notation SM corresponds to predictions of the 
Bahcall-Pinsonneault standard solar model (BPOO) of ref. [|l4|] and the standard model of 
electroweak interactions [^] (with no neutrino oscillations). The SNO rates (pure D2O 
phase) in the column labeled Result/SM are obtained from the published SNO spectral 
data by assuming that the shape of the 8 B neutrino spectrum is not affected by physics 
beyond the standard electroweak model. However, in our global analyses, we allow for 



spectral distortion. The SNO rates (salt phase) are not constrained to the B shape [33]. 



Experiment 



Observable (# Data) 



Measured/SM Reference 



Chlorine 

SAGE+GALLEX/GNOt 
Super-Kamiokande 
SNO (pure D20 phase) 



SNO (salt phase) 



KamLAND 



Average Rate (1) 
Average Rate (1) 
Zenith Spectrum (44) 
Day-night Spectrum (34) 



Average Rates (3) 



Spectrum (13) 



[CC]=0.34±0.03 
[CC]=0.54±0.03 
[ES]=0.465 ±0.015 
[CC]=0.35±0.02 
[ES]=0.47±0.05 
[NC]=1.01 ±0.13 
[CC]=0.32±0.02 
[ES]=0.44±0.06 
[NC]=1.03±0.09 
[CC]=0.61±0.09 



§ I M 



MM 



MM 



MM 
Hi 
Hi 

Hi 
PI 



t SAGE rate: 69.1±|;| ±|f SNU 



]; GALLEX/GNO rate: 69.3 ± 4.1 ± 3.6 SNU 



3.2 Global x 2 , free parameters, and reduced solar neutrino fluxes 

We calculate the global x 2 by fitting to all the available data, solar (83 measure- 
ments) plus reactor (13 measurements). Formally, the global \ 2 can be written in 
the form |34, |35 



Xglobal - Xsolar(A m > #12, {£, /b, /Be, fp-p, /cNo}) 

+ XKamLAND^m 2 , #12, #13) + XcHOOZ+ATM (#13) • (3-1) 

Depending upon the case we consider, there can be as many as nine free param- 
eters in Xsoiar> including, Am 2 ,6 12 , £ , / B , fs e , fp-p, and / C no (3 CNO fluxes, see 
below). The neutrino oscillation parameters Am 2 , 612 have their usual meaning. The 
reduced fluxes /b, /Be, fp-p, and /cno are denned below. The parameter £ describes 
departures from two-neutrino oscillations. In general, £ is a vector defined by several 
parameters. If we consider the standard scenario of three flavors, the higher mass 
squared averages out and £ = #13. 

The function XKamLAND depends only on Am 2 , 612 and #13. We discuss in section 
2.2 of ref. [0] our procedure for analyzing the data from the KamLAND experiment. 
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Figure 2: Ax 2 dependence on #13. Here #13 is the mixing angle between the first and 
the third neutrino eigenstates. The dependence of x 2 upon #13 that is shown in the figure 
corresponds to the global analysis of solar+atmospheric+KamLAND+K2K+CHOOZ data 
available in 2002 and displayed in Fig. 3c of ref. p3|] . We have used the 2002 dependence 
shown in the figure in calculating the global \ 2 defined by equation 3.1. We have not 
assumed a flat distribution for #13 below the CHOOZ bound. If a 1% p — p solar neutrino 
experiment is performed and analyzed in conjunction with a 5% 7 Be solar neutrino exper- 
iment, three years of KamLAND reactor data, and the currently available solar neutrino 
data, then the bound in the curve labeled ' + p — p ' can be achieved. 



For most of the calculations described in this paper, we marginalize Xgiobai making 
use of the function Xchooz+atm($i3) that was obtained in the updated analysis p4 



of atmospheric P5[ , K2K accelerator |57J, and CHOOZ reactor |J8], [39] data (see also, 
refs. [^2[ p3fl ). Figure |2| (which is the same as Fig. 3c of ref. IffiH ) shows the specific 
dependence upon 6*13 that we use. This formulation takes account of all the available 
data from solar, reactor, and accelerator neutrino oscillation experiments. We have 
not assumed, as is often done, a flat probability distribution for all values of #13 
below the CHOOZ bound. The fact that we take account of the actual experimental 
constraints on 6*13 decreases the estimated influence of 6*13 compared to what would 
have been obtained for a flat probability distribution. 

We have included in our analysis the possibility of an active-sterile admixture 
and in that case £ = 77, where 77 describes the possibility that v e oscillates into a state 
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that is a linear combination of active (i/ a ) and sterile (u s ) neutrino states, 



u e — > cos 77 z/ a + sin r\ v s 



(3.2) 



This admixture arises in the framework of 4-z/ mixing | 40| . The analysis procedure 
we follow in the case of a non-zero admixture of sterile neutrinos is described in 
section 3.3.2 of ref. 0. For this case only, we assume 9 X3 = 0. 

We follow the usual convention and represent the reduced 8 B solar neutrino flux 
by the parameter /b here /b is the ratio of the true flux to the flux predicted by the 
standard solar model. Thus 

0( 8 B) 



h = 

Here 0( 8 B)bpoo = 5.05 x 10 6 cm-V 1 
f Be , by analogy with f B . Thus 

fBe = 



s Pl • (3-3) 

O J BPOO 

We define the reduced 7 Be neutrino flux, 



7 Be) 



( 7 Be) 



BPOO 



(3.4) 



where 0( 7 Be)spoo = 4.77xl0 9 cm 2 s 1 with a la uncertainty of ±10%. The reduced 
p — p neutrino flux is defined by 



fp-p 



<P{P - P) 
<f)(p - p) 



BPOO 



(3.5) 



where |L4| 4>{p — p)bpoo = 5.95 x 10 10 cm _2 s~ 1 with a la uncertainty of ±1%. 

Similarly, we define reduced 13 N, 15 and 17 F neutrino fluxes that we parameterize 
by /cno- We have checked that different limiting cases in the relative weights of the 
13 N, 15 and 17 F fluxes (see discussion in ref. ||41|| ) do not affect the accuracy of our 
results. 



3.3 Survival probability and #13 

The survival probability of electron neutrinos in the case of three neutrino oscilla- 
tions, P e 3 e ", can be simply related to the survival probability, P^, for two neutrino 



oscillations by the equation [42 



cos 4 6 13 P^(Am 2 , 6>i 2 ; n e — > cos 2 9i 3 n e ) + sin 4 9 13 . 



(3.6) 



The effect of AM 2 , the mass difference squared characteristic of atmospheric neu- 
trinos, averages out in equation |3.6| for the energies and distances characteristic of 
solar neutrino propagation. The effective two-neutrino problem is solved with a 
re-normalized electron density cos 2 9\ 3 n e . The results from the CHOOZ reactor ex- 
periment pHl B9f place a strong upper bound on sin 2 2#i3 , implying that # 13 is close 



to or close to n/2. Atmospheric and solar data select the first option (cos 4 #13 
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close to 1 and sin 4 #13 close to 0). Thus the main effect of a small allowed #13 on the 
survival probability is the introduction of the factor cos 4 6*13 in equation [3~6. 



In what follows, we shall use 'survival probability' and P ee to denote P e 3 e ". Where 
numerical results are reported, we marginalize over 6^. At all points in oscillation 
parameter space, we use the value of 6*13 that minimizes x 2 for that set of parameters. 

3.4 Marginalization for predicted event rates 

In order to calculate properly the allowed ranges of the reduced event rates for 
[ 7 Be]„_ e (defined in equation |6.1| ) or [p — p] u - e (defined in equation [7.1|) , we first 



evaluate a global minimum \ 2 - We then calculate, for example, x 2 and [ 7 Be] v _ e 
at each of many points in the mult i- dimensional input parameter space. From this 
data set, we construct the minimum value of x 2 & t each value of [ 7 Be] !y _ e . Thus we 
are able to construct a function ;^ n ([ 7 Be]i,_ e ). We determine the allowed range of 
[ 7 Be] I/ _ e by selecting the values of [ 7 Be] y _ e that produce the appropriate Ax 2 for the 
specified confidence limit and 1 d.o.f. 

3.5 Simulated KamLAND data 

We simulate three years of KamLAND measurements using the data provided in 



the KamLAND paper ||18|| . We multiply by nine the number of events for all the 
prompt energy bins above 2.6 MeV. Following the KamLAND collaboration, we 
limit our analysis to energy bins above 2.6 MeV in order to eliminate the geo-neutrino 
background. The energy bins above 6 MeV contain zero events in the first KamLAND 
results. For these bins, we simulate the expected signal in 3 years of KamLAND data 
by using the number of events predicted by the present best fit oscillation solution 
(tan 2 6>i2 = 0.42, Am 2 = 7.1 x 10- 5 eV 2 ) for all the available solar plus KamLAND 
data. 

3.6 SNO salt phase data 



In the analyses presented in this paper, we have included the salt-phase SNO data |]33 
and the refined GNO g(J and SAGE § data, all released on September 7, 2003. Ev- 
ery one of the analyses performed in this paper was carried out with (in the current 
version of the paper) and without (in the previous version, hep-ph/0305159,v2) the 
new SNO, GNO, and SAGE data. Performing the analyses with and without the 
new data allowed us to double check all our conclusions and to verify that the depar- 
ture from maximal mixing was the only physical variable whose value or statistical 
significance was much affected by the addition of the new data 2 . 



2 In the version of this paper that was posted on the archive on June 19, 2003 and that was 
submitted for publication prior to the availability of the recently released SNO salt-phase data, we 
made best-estimate simulations of the expected salt-phase measurements (see hep-ph/0305159,v2, 



section 3.6). We shall see in section 3.2 that these simulations turned out to be accurate and the 
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In this version of our paper, we have included the charged current, neutral cur- 
rent, and elastic scattering rates |33| as determined by the SNO collaboration, where 
the spectral distributions of electron scattering and charged current events are not 
constrained by the known shape of the 8 B neutrino energy spectrum. We treated 
the errors in the salt phase measurements as described in the detailed user guide 
written by the collaboration [£3|; we no longer need to use our simulations of the 
SNO salt-phase data. 



4. Solar plus KamLAND constraints 

What can we say today about neutrino oscillation parameters and solar neutrino 
fluxes using a global analysis of all the existing data on solar and reactor neu- 
trino experiments? And, what improvements can we expect in our knowledge of 
the oscillation parameters and neutrino fluxes after three years of data taking by the 
KamLAND collaboration? 



Using a global analysis, we determine in section |4 . 1| the allowed regions of neu- 
trino oscillation parameters and the 8 B solar neutrino flux assuming that the 8 B 
neutrino flux is a free parameter. In this subsection, we assume that all solar neu- 
trino fluxes except the 8 B neutrino flux are described well by the standard solar 



model calculation [Tj] of their best-estimates and uncertainties. Allowing both the 
7 Be and the 8 B neutrino fluxes to be free parameters, we repeat the global analysis 
in section [4.2|. We discuss in section |0 the changes in the predictions for future 



solar neutrino experiments that are caused by allowing the 7 Be neutrino flux to vary 
freely. 



The most dramatic result is obtained in section [4.3| . In this subsection, we impose 
the luminosity constraint while allowing the p — p, 7 Be, 8 B, and CNO neutrino fluxes 
to vary freely. The existing solar neutrino and KamLAND experiments constrain 
the total p — p flux to ±2% (la)and are in agreement with the standard solar model 



prediction within the calculated error. 

We present in section |4.4| the constraints on the sterile neutrino fraction of the 
8 Be neutrino flux that are implied by the existing solar and KamLAND data. We 
also evaluate the effect of three years of simulated KamLAND data on the expected 
bound for the sterile fraction of the 8 B neutrino flux. 



conclusions that were deduced from the simulations are in agreement with those derived by the 
SNO collaboration from the actual measurements. The agreement between the simulations and 
the measurements for the SNO salt-phase provides circumstantial evidence that the simulations we 
make later in this paper for other experiments actually capture the essential features of those new 
experiments. Readers interested in simulating solar neutrino measurements may profit by reading 
section 3.6 of hep-ph/0305159,v2 to see what are the essential ingredients of an accurate simulation 
of the SNO data. 
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Figure 3: Allowed oscillation parameters: Solar plus KamLAND measure- 
ments. The figure shows the currently allowed la region for oscillation parameters 
that is obtained by a global fit, as described in the text, of all the available solar 
(S) H [?], ||, 15, [IB], [17], |2pJ, UK plus reactor (K) data. At 3a, the allowed regions form 
two connected islands |2|. The current best-fit point is labeled by a solid star. The figure 
also shows how the allowed contour is expected to be reduced after three years of operation 
of the KamLAND reactor experiment. In constructing these two contours, the 8 B solar 
neutrino flux was treated as a free parameter and all the other solar neutrino fluxes, and 
their theoretical uncertainties, were taken from the standard solar model of BP00 [14] The 
figure also shows a contour labeled 'B + lum', which corresponds to the case summarized 
in table [|, in which all the neutrino fluxes are treated as free parameters subject to the 
luminosity constraint. The la contours were drawn in all cases by marginalizing over #13 
and the 8 B solar neutrino flux. 



We investigate in section [4.5| the extent to which solar neutrino and reactor data 

we 



can constrain the CNO contribution to the solar luminosity and in section 4.6 



show that neutrino observations alone can provide a 20% measurement of the total 
solar luminosity. 

All our of investigations in this section show that we will be stuck with a large un- 
certainty in the experimentally-determined 7 Be solar neutrino flux until a dedicated 
7 Be solar neutrino experiment is performed. 

For the reader who is mainly interested in the bottom line, we recommend jump- 
ing directly to section fO. 
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Table 2: Allowed neutrino parameters and 8 B solar neutrino flux: Solar plus 
KamLAND measurements. For Am 2 , tan 2 #12, and the 8 B solar neutrino flux, the table 
presents the global best fit and la (3a) allowed ranges. The first row shows the allowed 
ranges determined from a global analysis of all existing solar neutrino (S')data available 
before September 7, 2003 as well as the initial KamLAND (K) data. The the second and 
third rows are determined from a global analysis of all existing solar (S) and KamLAND 
(K) data, including the SNO, GNO, and SAGE data @, ||, ||] released on September 7, 



2003. The 8 B solar neutrino flux is treated as a free parameter, but the 7 Be and all other 



solar neutrino fluxes are assumed to have the standard solar model 14] predicted values 
and uncertainties. 

4.1 8 B neutrino flux as a free parameter 

Figure |3] shows the la allowed region for the neutrino oscillation parameters Am 2 
and tan 2 # 12 that is obtained by a global fit to all of the existing solar plus KamLAND 
data. In addition, the figure shows the reduction in the allowed region that may be 
expected if the KamLAND experiment takes data for a total of three years. 

Table |2| shows the allowed ranges at la (3a) of the 8 B solar neutrino flux and 
the neutrino oscillation parameters Am 2 and tan 2 #12. In table |2|, we follow the usual 
convention and represent the 8 B solar neutrino flux by the parameter /b, where /b is 
the ratio of the true flux to the flux predicted by the standard solar model as defined 



in eq. [O 



The comparison of the first and second rows of table |] shows the impact of 



the recent SNO f33|, GNO [|(J, and SAGE data, especially the recently released 
SNO salt phase data |33|, |43[| . The positive 3a uncertainty in Am 2 is decreased by 



a factor of 4.5. The solar data before (after) the SNO salt phase release disfavors 
the higher LMA island less (more) than A% 2 = 9 . The SNO salt phase data also 
improved significantly our knowledge of tan 2 8u, especially the higher values; the 3a 
uncertainty in tan 2 (9 12 is decreased by a factor of order 1.5. Moreover, the la and 
3a uncertainties in /b are decreased by a 25%. 

The la ranges of Am 2 and tan 2 612 are slightly different in figure ^| and in table [| 
Figure ^ was constructed by marginalizing over # 13 and /b; the confidence level for 
the figure was calculated for two degrees of freedom. In constructing each column 
of table 0, we marginalized over all of the variables except the one whose range is 
listed; the confidence levels for the table were calculated for one degree of freedom. 

The longer operation of the KamLAND experiment will shrink significantly the 
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Table 3: Allowed neutrino parameters with free 7 Be and 8 B solar neutrino 
fluxes: Solar plus KamLAND measurements. The table presents the global best-fit 
and la (3a) allowed ranges for the 7 Be and 8 B solar neutrino fluxes and for the neutrino 
oscillation parameters Am 2 and tan 2 612- The allowed ranges are determined from a global 
analysis, with one degree of freedom, of all currently available solar (S) || ||, [H], 16, 17, 



2?J, H and KamLAND (K) |j data. The 7 Be and 8 B solar neutrino fluxes are treated 



as free parameters, but all other solar neutrino fluxes are assumed to have the standard 
solar model (BP00) predicted values and uncertainties. In constructing each column of the 
table, we have marginalized over all variables except the one whose range is shown. 



uncertainties in Am 2 . The la (3a) uncertainty in Am 2 is expected to decrease by a 
factor of order 2 (2.5). The expected improvements in our knowledge of tan 2 # 12 and 
of the 8 B neutrino flux are much more modest. 

4.2 7 Be and 8 B neutrino fluxes as free parameters 

Table ||] shows how the allowed range of the oscillation parameters and the solar 
neutrino fluxes change when the flux of 7 Be solar neutrinos is treated as a free 
parameter. Figure |5] shows the la contours in the allowed oscillation parameter 
plane for the case in which the 7 Be solar neutrino flux is treated as a free parameter 
(dashed curve) and also for the case in which the 7 Be neutrino flux is constrained to 
have the standard solar model best-fit value and uncertainty (dotted curve). These 
two extreme treatments of the 7 Be solar neutrino flux produce very similar allowed 
regions in the Am 2 -tan 2 # 12 plane. 

The allowed ranges of Am 2 , tan 2 #12, and /b are only marginally affected by 
letting the 7 Be neutrino flux vary. 

However, the first row of table |3] shows that the current la experimental uncer- 
tainty in the reduced 7 Be neutrino flux, /ee (see eq. |3.4| ), is a factor of 2.3 times larger 
than the quoted theoretical uncertainty in the standard solar model calculation of 
the 7 Be neutrino flux. Moreover, table [3] shows that the experimental uncertainty in 
determining the 7 Be solar neutrino flux is not expected to decrease significantly as a 
result of running the KamLAND experiment for three years. The KamLAND experi- 
ment does not provide a strong constraint on the oscillation probability for neutrinos 
with energies comparable to the 0.86 MeV possessed by the 7 Be solar neutrinos. 

We have carried out global solutions with and without including the measure- 
ment of the solar neutrino capture rate in chlorine (cf. table [5]). The purpose of 
these calculations was to determine how sensitive to a single experiment are the cur- 
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Table 4: Allowed neutrino parameters with free p— p, 7 Be, and 8 B solar neutrino 
fluxes: with and without luminosity constraint. Analysis A corresponds to a global 
analysis in which the p—p, 8 B, and 7 Be solar neutrino fluxes were treated as free parameters. 
For Analysis B, the CNO neutrino fluxes were also treated as free parameters as well as 
the p — p, 8 B, and 7 Be fluxes. The symbol '+ lum' indicates that for rows two and three 
of table H the luminosity constraint [|TJ] is included in the analysis. The table presents 
the global best fit and la (3a) allowed ranges for the p — p, 7 Be, and 8 B solar neutrino 
fluxes and for the neutrino mixing parameter tan 2 612- For all cases presented in this table, 
Am 2 = 7.3i°J x 10" 5 eV 2 The results given here were obtained using all the currently 
available data from the solar [§, 0> !> Hi © 0' M> IH and KamLAND ||] neutrino 



experiments. All other (much less important) solar neutrino fluxes are assumed to have 
the standard solar model (BP00) predicted values and uncertainties. In constructing each 
column of the table, we have marginalized over all variables except the one whose range is 
shown in the column of interest. 



rent inferences regarding neutrino oscillation parameters and solar neutrino fluxes. 
In our view, the implications of solar neutrino experiments are too important to 
rest solely on the results of a single measurement. Including the chlorine experi- 
ment causes only a modest improvement in the allowed ranges, of order 10%, for the 
neutrino oscillation parameters and for the 8 B neutrino flux. However, the chlorine 
measurement affects significantly the best-fit value for the 7 Be neutrino flux and 
the inferred uncertainty in this flux. With the chlorine experiment, f& e = 0.64^q;22 
(iaeit)) while without the chlorine measurement fs e = 0-97to'%2 (^92) ■ The chlorine 
measurement drives down the best-estimate for / Be by about 30% and decreases the 
la error on / Be by almost 30% . 



4.3 All neutrino fluxes as free parameters plus luminosity constraint 

Table |] presents the results of global analyses in which the p — p solar neutrino flux, 
as well as the 7 Be and 8 B solar neutrino fluxes, are treated as free parameters. This 
is the first table in which we show results obtained with and without imposing the 



luminosity constraint [12]. Cases in which the luminosity constraint are included are 



denoted by '+ lum' in the table. 



Figure f| illustrates (primarily) the impact of the SNO salt-phase data |33| 



an- 



nounced at the TAUP03 conference on September 7, 2003. The left panel of the 
figure shows the allowed oscillation regions (at 90%, 95%, 99%, and 99.73% (3a) ) 
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Figure 4: "Before and After TAUP03." The left panel shows the allowed regions 
for neutrino oscillations computed using all solar neutrino and reactor experimental data 
available prior to September 07, 2003. The right panel shows the allowed regions computed 
with the same procedures but including the SNO salt-phase data |33| and the improved 
SAGE and GALLEX-GNO @ data that were announced at the conference TAUP03 on 
September 7, 2003. The CL contours shown in the figure are 90%, 95%, 99%, and 99.73% 
(3a). The global best-fit points are marked by a star. 



before the salt-phase data were announced and the right panel shows the currently 
allowed oscillation regions. 

4.3.1 Neutrino fluxes without the luminosity constraint 

The first row of table which was calculated without including the luminosity 
constraint, should be compared with the second row of table |2| and the first row of 
table H The accuracy with which the global solution appears to determine the solar 
neutrino fluxes degrades as more free neutrino fluxes are varied in the analysis. For 
example, the quoted BP00 error on the predicted 7 Be neutrino flux is ±10% (lcr), 
whereas the experimental determination is only accurate to ±26% when the 7 Be 
neutrino flux is treated as a free parameter. Similarly, the SSM uncertainty on the 
predicted p — p neutrino flux is ±1%, but the flux is only determined experimentally 
to ±21%. The 8 B solar neutrino flux, which is primarily determined by the SNO 
and Super-Kamiokande measurements, is not affected much by treating additional 
neutrino fluxes as free parameters. The accuracy with which Am 2 and tan 2 812 are 
determined is also not affected significantly by adding more free fluxes. 

The luminosity constraint is grossly violated if we allow the p — p, 7 Be, and 8 B 
neutrino fluxes to vary as free parameters without requiring that the total luminosity 
equal the observed luminosity. The first row of table [| shows that the best-fit value 
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for the p — p neutrino flux is / p _ p = 1.38 in this case, whereas the maximum value 
allowed by the luminosity constraint is / p _ p = 1.09 (see equation |5~l"l) . 



4.3.2 Neutrino fluxes with the luminosity constraint 

The constraints on the p — p neutrino flux are dramatically improved, by more than 
a factor of ten, when the luminosity constraint is applied. This improvement is 
independent of whether the CNO neutrino fluxes are assumed to have their standard 
solar model best-fit and uncertainties (row two of table [|) or whether the CNO fluxes 
are treated as free parameters (row three of table f§). Remarkably, the p — p solar 
neutrino flux is known at present to ±2% (la uncertainty) and is in agreement to 
that accuracy with the prediction of the standard solar model 3 

The constraints on the 7 Be neutrino flux are not improved by imposing the 
luminosity constraint. The luminosity constraint affects differently the p — p and 
7 Be solar neutrino fluxes. We shall explain the origin of this difference in section |[ 
However, the alert reader may notice a counter-intuitive result shown in column 4 
(/Be) of table f|. The uncertainties are actually larger for f Be when the luminosity 
constraint is imposed. The reason for this anomalous behavior is that the global Xmm 
is smaller (Xmm = 78.6) for case A of table f| (without the luminosity constraint) 
than the global x min (80.4) found for case A plus the luminosity constraint. The 
uncertainties at a fixed confidence level (Ax 2 ) are calculated by computing the range 
of f Be allowed by the relation Ax 2 = X 2 (/b c ) - Xmm- 

The CNO fluxes are poorly constrained by the available solar neutrino data (cf. 
ref. [f|I[]). In all cases, that we have considered the favored solution has the 13 N, 15 O, 
and 17 F solar neutrino fluxes all equal to zero. The 15 flux is the best determined 
of the CNO solar neutrino fluxes. We find in units of the standard solar model 



prediction [J E that 



fo = 0.0^0. (4.1) 

The constraint given in equation JO] is unaffected, to the accuracy shown, by whether 
or not the luminosity constraint is imposed. 

At la, the 13 N solar neutrino flux could be as large as 13 times the standard 
solar model prediction and the 17 F neutrino flux could be as large as 220 times the 
solar model prediction. 

3 The results given here should be compared with an insightful analysis given in section VI of 
ref. 0. The authors of ref. made a series of simplifying assumptions in analyzing the pre- 
KamLAND solar neutrino data and obtained / P - P = 1.29 ± 0.30. The principal assumptions 
made in this analysis, which did not include a global treatment of all the experiments together, 
included assuming the correctness of the BPOO fluxes for the 7 Be, CNO, and pep neutrinos and the 
constancy of the survival probability for all of these neutrinos. Although the authors suggest that 
their uncertainty may be underestimated, their reasoning shows clearly how the gallium experiments 
constrain the p ~ p neutrino flux. 
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4.3.3 Allowed oscillation solutions: with and without the luminosity con- 
straint 

In ref. @, the authors found that, using the existing solar plus KamLAND data, only 
the LMA oscillation solution is allowed at 4.7 a. Repeating the same analysis but 
including the SNO salt phase data and the other new solar neutrino data 0, |33|, pOj ], 
we find that only the LMA oscillation solution is allowed at 5.2 a. 

All other solar neutrino oscillation solutions (e.g., LOW, VAC, SMA) are ex- 
cluded. In the analysis described in ref. and in other similar discussions (see, e.g., 
ref. ||), the solar neutrino fluxes and their uncertainties were either all taken from 
the standard solar model or all of the fluxes except the 8 B neutrino flux were taken 
from the standard solar model. 

If, instead, all solar neutrino fluxes are regarded as free parameters and the 
luminosity constraint is imposed, we find that the LMA oscillation is the only allowed 
solution at 4.5 a. 

All solutions in which the neutrino energy spectra are undistorted, i.e., all so- 
lutions without some form of neutrino oscillations (what used to be called 'non- 
standard solar models'), are disfavored at 8.3a fll2 |. 



The existing constraints on the 7 Be and p — p neutrino fluxes result from the 
chlorine || and gallium || ^Of solar neutrino experiments (cf. table |5|). It is of 



interest to test separately the sensitivity of the inferred constraints on the neutrino 
fluxes by omitting either the chlorine or the gallium experiments from the global 
analysis (but imposing the luminosity constraint as in the third row of table 
Omitting the chlorine experiment, we find that the constraints on Am 2 , tan 2 812, 
and /b are essentially as stated in the third row of table ^, but that the constraints 
on the 7 Be and p — p neutrino fluxes are much relaxed: 

Without CI: / Be = 1.24±?;|1 ; f p . p = l.OOi^CT- (4-2) 

Omitting the gallium experiment, we find 

Without Ga: / Be = 0.33l*J* ; f p . p = 1.05±™ 2 (^)- (4-3) 

The chlorine experiment has a stronger effect on the constraints than does the 
gallium experiment. 

Figure |3] shows that the uncertainties in determining Am 2 and tan 2 812 from 
existing data are not affected significantly by whether only the 8 B neutrino flux 
is treated as a free variable or whether, instead, all the fluxes are treated as free 
variables subject to the luminosity constraint. In both cases, the allowed ranges are 
practically the same (cf. the second row of table |2] with the last row of table f|). 

4.4 Sterile neutrinos 

We can parameterize the sterile contribution to the neutrino flux in terms of sin 2 r\ 
or, alternatively, in terms of a derived parameter /b sterile, the sterile fraction of the 
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8 B neutrino flux (see discussion in section ^2 and the more complete description in 
ref §). 

The la allowed range for the active-sterile admixture is 



sin 2 77 < 0.09 (4.4) 

in our analysis of the existing solar plus KamLAND data (cf. the first row of table |2|). 
The fundamental parameter describing the sterile fraction is sin 2 r\ (cf. section [O] 
and ref 0). However, it is convenient to think in terms of a sterile fraction of the 
flux, /b, sterile , which is potentially observable in the Super-Kamiokande and SNO 
experiments. This range corresponds to 

/b, sterile = O.Ot°;°o (solar + KamLAND) . (4.5) 

Three years of accumulation of data by KamLAND (the S + K 3 yr second row of 
table 0) can marginally improve this bound to 

sin 2 7] < 0.08, (4.6) 



or 

/ B ,steriie = 0.0l°;°(j (solar + KamLAND 3 yr) . (4.7) 

On the other hand, if we impose the luminosity constraint while allowing the 
p — p, 7 Be, 8 B, and CNO neutrino fluxes to vary freely (analysis B + lum in table ^), 
our present bounds are slightly weaker. In this case, we find 

sin 2 7/<0.10, (4.8) 



or 



/b, sterile = 0-Olaoo (solar + KamLAND; free fluxes; luminosity constraint) . (4.9) 
4.5 CNO luminosity 

How much of the solar luminosity is due to the carbon-nitrogen-oxygen (CNO) nu- 
clear fusion reactions that Hans Bethe suggested in 1939 are primarily responsible 
for energy generation in the Sun? 



In ref. |)4T |, the authors derived la (3a) constraints on the allowed range of 
the CNO contribution to the solar luminosity using all the presently existing solar 
neutrino and KamLAND experimental data. Imposing the luminosity constraint and 
carrying out a global analysis of the ratio of the CNO-produced luminosity to the 
total luminosity, they found pH| 

= o.Oijj-{j % (O.Oj£o %) at la (3d) (with luminosity constraint) . (4.10) 
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We have carried out a similar global analysis but using simulated three years of data 



for KamLAND. The results summarized in equation |4.10| are not affected, to the 
accuracy shown, by the additional simulated KamLAND data. 

For completeness, we have also carried out a global analysis allowing the fluxes 
to be free but without imposing the luminosity constraint. We find 

= 1.6t\'t% (l-6lio%) at 1(7 (3a) (no luminosity constraint) . (4.11) 

Lq 

The best-estimate values given in equation |4.10| and equation |4.11| are different. 



If the luminosity constraint is not imposed, then the global solution prefers a non- 
zero value for the CNO luminosity that is just slightly larger than the standard solar 
model best-estimate of 1.5%. The low-energy neutrinos that are required to fit the 
radiochemical experimental data are provided in this case by the CNO reactions. 
Moreover, the chlorine and gallium experiments suggest a preferred value of zero for 
the 7 Be solar neutrino flux (see case A of table §). 

However, if the luminosity constraint is imposed, then the 7 Be solar neutrino 
flux cannot easily be zero and the required low energy neutrinos are produced by 
7 Be electron capture (see case 'B + lum' of table HI and the approximate form of 



the luminosity constraint, equation STT). In this case, the radiochemical experiments 
prefer a zero-value for the CNO luminosity. 

The upper limits given in equations |4.10| and [4.11| are relatively insensitive to 



whether or not the luminosity constraint is imposed. 
4.6 Total luminosity measurement with neutrinos 

What limits do existing solar plus reactor neutrino experiments place on the total 
luminosity of the Sun? 

Each solar neutrino flux, $j, corresponds to a well-defined flux of energy, a^i, 



where the energy coefficients are given in table 1 of ref. [12|]. This energy is, of 
course, much larger than the actual energy of the neutrinos since most of the energy 
from solar fusion reactions is released to the star in the form of thermal energy and 
eventually radiated away in the form of photons (see discussion in ref. [|T2|). 

We have performed a global analysis of all the available solar and reactor data to 
find the la (3a) allowed range for the total flux of solar luminosity, based 
solely upon the neutrino data. We find for the ratio of the neutrino-inferred solar 
luminosity, L & (neutrino — inferred), to the accurately measured photon luminosity, 
L Q , that 

L Q (neutrino — inferred) r +0 2 , +0 7 



1-4^ (!"oI ■ (4-12) 

At 3cr, the neutrino-inferred solar luminosity can be as large ctS (ctS small as) 2.1 
(0.8) the precisely measured photon-luminosity. Three years of KamLAND data is 
expected to improve the limits stated in equation |4.12| by only about 10%. 
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Table 5: Neutrino oscillation predictions for the chlorine and gallium radio- 
chemical experiments. The rates without parentheses (in parentheses) are presented 
for the best-fit oscillation parameters of the allowed solutions listed in table |2| , which as- 
sumes neutrino fluxes from the BPOO solar model except for the 8 B neutrino flux (table 
which was derived with free neutrino fluxes and the luminosity constraint). The predictions 
are based upon the global analysis of existing solar plus KamLAND data and the neutrino 
absorption cross sections given in refs. |p9| , |44| , jig ]. The total rates should be compared 
with the standard solar model values flTJ], which are 7.6ii;f SNU (chlorine) and 1281? 
SNU (gallium), and the measured values, which are 2.56 ±0.23 (chlorine ||) and 69.2 ±4.0 
(gallium @, !,§]). 



Source 


CI 


Ga 




(SNU) 


(SNU) 




LMA 


LMA 


p — p 


0(0) 


40.6 (41.8) 


pep 


0.13(0.13) 


1.57 (1.61) 


hep 


0.02(0.02) 


0.03 (0.03) 


7 Be 


0.64(0.60) 


19.0 (17.9) 


8 B 


2.02(2.00) 


4.28 (4.23) 


13 N 


0.05(0.00) 


1.83 (0.00) 


15Q 


0.18(0.00) 


2.94 (0.00) 


17p 


0.00(0.00) 


0.03 (0.00) 



Total 3.04 ±0.04(2.75 ±0.18) 70.3±^(65.6j:£i) 



4.7 Predictions for the radiochemical chlorine and gallium experiments 

Table |5] summarizes the predictions for the radiochemical solar neutrino experiments 
of the LMA solution that are inferred from the existing solar neutrino and KamLAND 
experiments. The predictions are based upon the global analyses for which the best- 
fit neutrino oscillation parameters are listed in the first row of table |2| (and the third 
row of table f|, free solar neutrino fluxes plus luminosity constraint). The measured 
rate in the chlorine experiment is approximately 2cr (la) smaller than the best-fit 
LMA prediction. 

5. How does the luminosity constraint affect the allowed range 
of solar neutrino fluxes? 

The luminosity constraint states that a specific linear combination of all the neutrino 
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fluxes that are produced by fusion reactions at a temperature of order 1 keV must 
be equal to unity. The precise coefficients for this linear combination are given in 
table 1 and equation 2 of ref. [[12] . 

We can understand the strong effect of the luminosity constraint on the p — p 
neutrino flux, and the relatively weak effect on the 7 Be neutrino flux, by including 
only the largest terms in the luminosity constraint. We define a reduced neutrino flux 
to be the ratio of the true solar neutrino flux from a particular nuclear reaction to 
the neutrino flux predicted for this source by the BPOO solar model. If the standard 
solar model is a relatively good approximation to the actual Sun, then the leading 
terms in the luminosity constraint are 

(f){p-p) « 1.09 - O.O80( 7 Be) - O.O10( 15 O). (5.1) 



We see immediately from equation |5J] that the luminosity constraint does not allow 
the p — p flux to exceed, at any confidence level, 9% of the current standard solar 
model predicted flux. By contrast, if we do not impose the luminosity constraint, 
the best-fit value for <p(p — p) is 1.49 not 1.0 and the maximum value of <fi{p — p) is 
2.02 at 3<t (see the last column of tabled). The strong lower limit on <f)(p—p) follows 



from equation |5.1| because the chlorine and gallium solar neutrino experiments imply 
that 0( 7 Be) can not be much larger than unity (cf. table |5|). 

The small coefficient of 0( 7 Be) in equation |5.1| is the reason that the luminosity 
constraint does not impose tight constraints on the allowed 7 Be solar neutrino flux, 
provided that the standard solar model is a reasonable approximation to the actual 
Sun. In principle, the 7 Be solar neutrino flux could be as large as three times the 
predicted BPOO 7 Be neutrino flux without violating nuclear physics or energy con- 



straints |I2| . However, this would be true only for a real Sun that is very different 
from the standard solar model. 



6. What will we learn from a 7 Be solar neutrino experiment? 

In this section, we investigate quantitatively what we may expect to learn from a 
measurement of the 7 Be solar neutrino flux. 

We begin in section |6.1| by comparing the differences in the predictions for fu- 
ture experiments caused by treating the 7 Be and other solar neutrino fluxes as free 
parameters rather than by taking the fluxes, and their associated uncertainties, from 
the standard solar model. Then in section |6.2| we investigate how accurately a mea- 
surement of the 7 Be solar neutrino flux is expected to determine neutrino oscillation 
parameters. We study in section how accurately a 7 Be experiment will determine 
solar neutrino fluxes. We determine in section 16.41 how well a 7 Be neutrino-electron 



scattering experiment can constrain the CNO-produced luminosity of the Sun and 
in section |6.5| we show how accurately a 7 Be neutrino-electron scattering experiment 
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will help measure the total luminosity of the Sun. We compare in section 6J; what 
can be learned from a 7 Be CC experiment with what can be learned from a v — e 
scattering experiment. 

In all cases, we treat the 8 B solar neutrino flux as a free parameter. 

6.1 7 Be free: Implications for future experiments 

Whether or not / Be is treated as a free parameter affects strongly the uncertainty, 
and the best estimate, of the predicted rate at which electrons scatter 7 Be neutrinos 
or are absorbed in a CC experiment. We define the reduced 7 Be v — e scattering or 
CC (absorption) rate as follows: 

7 _ Observed rate , . 

^ ^ = BPOO predicted rate ' ^ ' ' 

where the denominator of equation ^71] is the rate calculated with the BPOO 7 Be 
neutrino flux assuming no neutrino oscillations. 

If we perform the global analysis of solar plus reactor experiments assuming 
that the BPOO calculation for the 7 Be solar neutrino flux and its uncertainty are 
valid, then the predicted rate in a 7 Be v — e scattering experiment is, with la (3<r) 
uncertainties: 



7 Be 



0.66 ± 0.02 (Zo l) ■ (6.2) 



Thus, if we use the solar model calculation of the 7 Be neutrino flux and its uncer- 
tainty, then the predicted event rate for the 7 Be rate experiment has a precision of 
±3%. The corresponding prediction for the reduced event rate in a CC (absorption) 
experiment with 7 Be solar neutrinos is: 



cc = 0.57 ± 0.02 (6.3) 



The precision of the CC prediction is comparable to what is obtained for a v — e 
scattering experiment. The results given in equation and equation pTSIonly make 



use of the currently available KamLAND data, but the results are practically the 
same if we use the three year simulated KamLAND result. 

The result shown in equation |3.3| for the CC reaction includes only 7 Be neutrinos 
from the 0.86 MeV electron capture line. On the other hand, the result for v — e 
scattering that is given in equation |6.2j includes neutrinos from all the solar neutrino 
fluxes that produce recoil electrons with energies in the range 0.25-0.8 MeV. Most 
of the recoil electrons in the selected energy range are produced by 7 Be neutrinos, 
although there are small contributions from CNO, pep, and p — p neutrinos (whose 
fluxes are taken from the BPOO solar model fUjl). 
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The situation is very different if we treat the 7 Be solar neutrino flux as a free 
variable in the global analysis. In this case, using the existing KamLAND data and 
all the available solar data, we find: 



r Be 



0.47 



+0.13/+0.39\ 
-0.12l-0.41 J 



and 



[free 7 Be flux] 



(6.4) 



r Be 



cc 



3g+ 0A4 f+ 0A1 ^ \f 7 



J -0. 131-0.38 J 



[free Be flux] . 



(6.5) 



The reason that the 3a lower limit in equation |6~4] and equation |0| does not reach 
zero is that the CNO contribution between 0.2 MeV and 0.8 MeV is included, by 
definition, in [ 7 Be]. 

The most important case allows all the neutrino fluxes to vary freely within the 
limits imposed by the luminosity constraint. With this experimental approach (cf. 
the third row of table f|), we find an even larger allowed range for the predicted rate 
of a 7 Be solar neutrino experiment. For this 'all free' case, we obtain 



7 Be 



0-58+^(^) [all free] 



(6.6) 



7 Be 



cc 



0.54+S-" 



h0.41> 
0.301-0.54; 



[all free] . 



(6.7) 



The uncertainties are enormous for the predicted rates in 7 Be experiments (cf. 
equations |6.6| - |6.7p if the solar neutrino fluxes used in the global fits to experiments 
are unconstrained by solar model calculations. The uncertainties in the predictions, 
~ ±60% of the best-estimate value for the u — e scattering rate (~ ±70% for the 
CC rate), are more than an order of magnitude larger than the uncertainties in the 
predictions, ~ ±3% for the u — e scattering rate (~ ±5% for the CC rate), if the 7 Be 
neutrino flux is constrained by the standard solar model prediction and uncertainties. 

6.2 How accurately will a 7 Be solar neutrino experiment determine neu- 
trino oscillation parameters? 

We suppose in this subsection that the u — e scattering rate, [ 7 Be]„_ e , is measured 
and that the best-fit value coincides with the result obtained assuming the standard 
solar model [14| 7 Be neutrino flux and the preferred values for Am 2 and tan 2 9\2 (cf. 



table g and equation |6.2|) . We have then carried out a global solution including all the 
currently available solar neutrino data, the simulated KamLAND three year data, 
and the hypothetical measurement of [ 7 Be] y _ e . We use the simulated three year 
data rather than the currently available one-year KamLAND data since it seems 
likely that the KamLAND reactor experiment will be completed before a 7 Be solar 
neutrino experiment is completed. 
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Experiments tan 2 9i 2 fse /i 



/'/' 





s + 


K 3 yr 


u - 4i -0.04 


/+0.15\ 
1-0.09J 


0.931 


-0.25 
0.63 


/+0.80\ 
1-0.93J 


1.02: 


-0.02 
-0.02 


'+0.06\ 
v-0.06j 


+ 


[ 7 Be], 


— e 


±30% 


U-^ 1 - 0.04 


1 +U.10N 

1-0.09^ 


0.991 


n on 

U.ZU 

0.26 


/_i_n t\~\ \ 
1-0.99J 


1.0151 


-U.Ulo 
-0.015 


/_i_n n/iv \ 

I-0.047J 


+ 


[ 7 Be] u 


— e 


±10% 


41 +0 - 04 

u -^ 1 - 0.04 


/+0.15N 
1-0.09^ 


1.071 


-0.10 
0.11 


/+0.31\ 
\-0.4b) 


1.0091 


-0.008 
-0.008 


(+0.024\ 
1-0.028J 


+ 


[ 7 Be} u 


— e 


±5% 


41 + - 04 

U-^ 1 - 0.04 


f+0.15\ 
1-0.09J 


1.091 


-0.05 
-0.07 


/+0.17\ 
1-0.31J 


1.0081 


-0.005 
-0.005 


f+O.Om 
1-0.024J 


+ 


[ 7 Be} u 


— e 


±3% 


41 + - 04 

U-^ 1 - 0.04 


/+0.15\ 
1-0.09^ 


1.091 


-0.04 
0.05 


/+0.12\ 
\-0.2Q) 


1.0081 


-0.002 
-0.005 


f +0.008 \ 
1-0.024J 


+ 


[ 7 Be}„ 


— e 


±1% 


41 + - 04 

U-^ 1 - 0.04 


f+0.15\ 
1-0.09J 


LlOj 


-0.02 
-0.05 


(+om\ 
1-0.24J 


1.0071 


-0.002 
-0.004 


/+0.005\ 
1-0.022J 



Table 6: The effect of a measurement of the 7 Be v — e scattering rate. The table 
presents best-fit and la (3<r) ranges computed by including all the currently available solar 
neutrino data (S) g 0, |, [TJ, [JJ|, [TJ, || plus three years of simulated KamLAND 



measurements subject to the luminosity constraint. The p — p, 7 Be, 8 B, and CNO neutrino 
fluxes are treated as free parameters. The last five rows of the table were computed using 
also a simulated measurement of the reduced v — e scattering rate (cf. equation |Q| ), 
[ 7 Be] i/ _ o = 0.66, assumed to be measured with 30%, 10%, 5%, 3% or 1 % precision. In 
constructing each column of the table, we have marginalized over all variables except the 
one whose range is shown. The allowed ranges of Am 2 and /b are, for all cases shown 
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in table §, Am 2 = 7.3±g | x 10~ 5 eV 2 and / B = 1.01±g;gj. The allowed range of tan 2 6 
changes only by a very small amount as a result of a precise 7 Be solar neutrino experiment. 



We have imposed the luminosity constraint \Y1\ on the solar neutrino fluxes 
when carrying out the global solution. We have treated the p — p, Be, 8 B, and CNO 
neutrino fluxes as free parameters. 

Table |6] compares the globally allowed ranges of tan 2 6*12, and fse if a 7 Be exper- 
iment is carried out with an accuracy of 30%, 10%, 5%, 3%, or 1% . For calibration, 
we also present in the first row of table |]the allowed parameter ranges in the absence 
of a 7 Be experiment. 

The allowed ranges of Am 2 and /b are, for all cases shown in table |6], Am 2 = 
7.3±° i x 10~ 5 eV 2 and f B = 1.01±g;g|. The allowed range of tan 2 12 changes only by 
a very small amount as a result of a precise 7 Be solar neutrino experiment. 

The solid curve in figure [5] represents the expected la allowed region after a 5% 
measurement of [ 7 Be] l/ _ e . This 'post- 7 Be' contour is not significantly smaller than 
the two other contours in figure ||]. The two other contours represent the result of 
analyzing with different options all the current solar neutrino data plus three years 
of simulated KamLAND data. The two options are to allow all the solar neutrino 
fluxes to be free parameters, subject to the luminosity constraint, or to treat only 
the 8 B solar neutrino flux clS db free parameter. 

We have also explored the effect of a 1% measurement of the 7 Be neutrino- 
electron scattering rate on the existing bound for the sterile fraction, sin 2 rj. We 
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- S + K 3 yr, 8 B free 

S + K 3 yr, all free 

+ [ ? Be] v _ e ± 5%, all free 




0.3 



0.4 



0.5 



0.6 



tan 2 
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Figure 5: Effect of 7 Be constraints on the allowed regions for oscillation param- 
eters. The main result illustrated in this figure is that the allowed region for oscillation 
parameters is not expected to be affected much by a measurement of the 7 Be solar neutrino 
flux. Each la contour in the plane of allowed oscillation parameters is constructed using 
all the currently available solar neutrino data |], |7], |8|, [l5|, [16|, [17], g0|, |3^] plus a simulated 
three years of KamL AND data (cf . figure |3|) , as well as a particular treatment of the solar 
neutrino fluxes. The dashed contour that extends to the smallest values of tan 2 9\2 corre- 
sponds to treating the 7 Be and all other solar neutrino fluxes as free parameters, while the 
contour that extends to the largest values of tan 2 #12 was computed by assuming that all 
of the solar neutrino fluxes except the 8 B neutrino flux are constrained by the estimated 
values and uncertainties obtained from the standard solar model 14 j. A 5% measurement 
of the v — e scattering rate for 7 Be solar neutrinos produces the solid intermediate curve 
in figure [5]. The contours were calculated by marginalizing over #13, /b, and /Be- 



find that even a 1% 7 Be experiment does not improve the existing limits derived in 



section 4.4 



We conclude that a measurement of the 7 Be solar neutrino flux will not contribute 
much to the knowledge of the allowed region in oscillation parameter space unless 
there is additional new physics beyond neutrino oscillations or unless one of the 
previous solar neutrino experiments has a large previously unrecognized systematic 
uncertainty or bias. 
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6.3 How accurately will a 7 Be neutrino experiment determine solar neu- 
trino fluxes? 



In the following subsection, section |6.3.1|, we study what can be learned about the 



7 Be solar neutrino flux by measuring the 7 Be v — e scattering rate. In section |6.3.2 



we investigate how much such 7 Be experiments will tell us about the p — p solar 
neutrino flux. 

6.3.1 What do we learn about the 7 Be solar neutrino flux? 

If [ 7 Be]j,_ c can be measured to an accuracy of ±10%, then, according to table || 
the experimental knowledge of the 7 Be solar neutrino flux will be improved by more 
than a factor of four over what will be known from the solar neutrino experiments 
plus three years of KamLAND reactor measurements. The experimental uncertainty 
in determining the 7 Be solar neutrino flux would then be essentially equal to the 
current uncertainty in the standard solar model prediction of the 7 Be neutrino flux. 
A 3% measurement of [ 7 Be] !/ c would lead to an improvement of a factor of eight in 
the accuracy with which the 7 Be solar neutrino flux is known. 

We conclude from a study of table ^ that a measurement of the 7 Be v — e 
scattering rate will lead to an enormous improvement in the experimental knowledge 
of the 7 Be solar neutrino flux. Without a specific 7 Be solar neutrino experiment, 
our empirical knowledge of the 7 Be solar neutrino flux will be a factor of four less 
precise than the current uncertainty in the standard solar model prediction [±40% 
(experimental) versus ±10% (standard solar model) ]. 

How sensitive are our quantitative conclusions to the assumed value for [ 7 Be]^_ e ? 
We have answered this question by carrying out global fits in which the assumed 
value for [ 7 Be] l/ _ e differs by ±lcr from the current best fit. The results obtained for 
the ranges of the neutrino oscillation parameters and the solar neutrino fluxes by 
assuming the la-different values are almost identical to the results shown in table |6] 
using our current best-prediction for [ 7 Be]^_ c . 

6.3.2 What do we learn about the p — p solar neutrino flux? 

We have already seen from the discussion of table f|that the p — p solar neutrino flux 
is already known to an accuracy of ±2% when the luminosity constraint is imposed. 
By comparing the first rows of table f| and table |6] that additional measurements 
by KamLAND are not likely to reduce greatly the uncertainty in the experimental 
measurement of the p — p flux. 

However, a 7 Be measurement can enormously improve our empirical knowledge 
of the p — p neutrino flux. A measurement of [ 7 Be] jy _ e to an accuracy of 5% will 
improve our experimental knowledge of the p — p flux by a factor of four, to an 
extraordinary ±0.5%. A 3% would reduce the experimental uncertainty to ±0.4% , 
a factor of five improvement. No significant further improvement would be achieved 
by a 1% measurement. 
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We investigate in section ^O] what constraints can be placed on the CNO- 
generated luminosity of the Sun, and in section |6.5| on the total solar luminosity, 
from a 7 Be neutrino-electron scattering experiment. 

6.3.3 What do we learn about the CNO solar neutrino fluxes? 

A measurement of the 7 Be solar neutrino flux would significantly improve our current 
very poor experimental knowledge of the CNO neutrino fluxes (cf. equation |4.1|) . 
Assuming a measurement of the 7 Be neutrino flux accurate to 5% and three more 
years of operation of KamLAND, we find that the 15 flux (in units of the standard 
solar model prediction [111]) could be constrained to be 



fo = 0.0±g^(±g:g) . (6.8) 

Thus a ±5% measurement of the 7 Be solar neutrino flux could improve the current 
la (3a) constraint on the 15 solar neutrino flux by a factor of five (two) (cf. equa- 
tion [4.1|) . At la, the 13 N solar neutrino flux could be as large as 6 times the standard 
solar model prediction and the 17 F neutrino flux could be as large as 40 times the 
solar model prediction. 

A measurement of the p — p solar neutrino flux would not by itself significantly 
improve our knowledge of the CNO fluxes. For example, adding the data from a 
simulated ±5% measurement of the p — p neutrino flux would improve the la upper 



limit given in equation 6.S from 0.5 to 0.4 and the 3a limit from 3.0 to 2.8 . 



6.4 What do we learn about the CNO luminosity? 

A measurement of the 7 Be neutrino-electron scattering rate to an accuracy of 5% 
(la) will constrain the CNO luminosity to 

LcNO -= O.OtH % (0.0±g;g %) at la (3a) (with 5% 7 Be measurement) . (6.9) 



The result shown in equation |6.9| is almost a factor of two improvement over the 
constraint that is obtained without a direct 7 Be measurement (cf. equation 4.10] 



and equation |4.11| ). The constraint on the CNO luminosity will not be significantly 
improved by making a 7 Be measurement that is more accurate than 5%. 

Whether or not the luminosity constraint is included does not affect, to the 



accuracy shown, the result given in equation 675 



6.5 How well can we constraint the total solar luminosity with a 7 Be 
neutrino experiment? 



Following the reasoning outlined in section |4.6| , we can calculate the constraints on 
the total from neutrino experiments, assuming a measurement accurate to 5% of the 
7 Be neutrino-electron scattering rate. We find 
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L Q (neutrino — inferred) 1 nf7+0 ,i3 f+o 4o\ reo/ 7 



1.071° 13 (±639) (5% 7 Be experiment) . (6.10) 
The result given in equation |6.10| is not changed, to the accuracy shown, if we assume 



either a 3% or a 1% measurement of the 7 Be neutrino-electron scattering rate. 

The 7 Be measurement as represented in equation |6.10| would improve by about 
50% the determination with neutrinos of the total solar luminosity (cf. equation [4.12|) . 

6.6 Will we learn more or less from a CC 7 Be experiment? 

We have concentrated in this section on v — e scattering experiments since the 



BOREXINO |46] and KamLAND [18[] detectors are already well advanced toward 
measuring the 7 Be neutrino-electron scattering rate. Moreover, there may be signif- 
icant uncertainties in the CC cross section that are not present for the v e scattering 
measurement. 

However, we note that a measurement with a CC detector with the same total 
experimental precision as for a v — e scattering experiment would give comparable, 
actually somewhat better, precision. For example, for a common measuring error of a 
measurement of ±5% (lcr), the 7 Be CC absorption rate would yield a determination 
of the total 7 Be flux with an error (±4.5%, la) that is about three-quarters of the 
error (±6%) on the total 7 Be neutrino flux that can be inferred from a v— e scattering 
experiment. In both cases, the global analysis includes a fit to all the available solar 
neutrino plus KamLAND data. 

7. What will we learn from a p — p solar neutrino experiment? 



To the best of our knowledge, Nakahata [47]] first discussed in a systematic way 
the possibility of using a measurement of the v — e scattering rate of p — p solar 
neutrinos to determine more accurately the mixing angle represented by tan 2 # 12 . In 
his discussion, Nakahata assumed that the total p — p neutrino flux is given by the 
standard solar model with its associated uncertainties (±1% for the predicted p — p 
neutrino flux) and therefore he did not apply the luminosity constraint. 

In what follows, it is convenient to use the reduced v — e scattering rate [or CC 
(absorption) rate] for p — p neutrinos which is defined by the relation: 

_ Observed rate (7 1) 

^ ~ ^ = BP00 predicted rate ' ^ ' ' 



where the denominator of equation 64. is the rate calculated with the BP00 p — p 
neutrino flux assuming no neutrino oscillations. 

We assume throughout this section that the 7 Be neutrino scattering rate has 
been measured to an accuracy of ±5% and that the data are available for three years 
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of operation of the KamLAND reactor experiment. From our current perspective, it 
seems likely that both the KamLAND reactor data and the 7 Be solar neutrino data 
will be available before the completion of a p — p solar neutrino experiment. 



In section we assume that the total p — p neutrino flux, and the uncertainty 
in the flux, are calculated accurately using the standard solar model. Following 
Nakahata, we explore what can be learned from a precise measurement of the rate, 
[p — p\ v -ei of neutrino-electron scattering of p — p solar neutrinos if we have faith in 
the standard solar model predictions. 

We consider a more rigorous and informative treatment in section |7.2| . In this 
analysis, the p — p, 7 Be, 8 B, and CNO neutrinos are treated as free parameters in a 
global analysis of all the data, with the data including simulated measurements of the 
p — p and 7 Be neutrino-electron scattering rates. We apply the luminosity constraint 
in the free-flux analysis described in section |7.2| , obtaining simulated powerful limits 
on the p — p neutrino flux and the range of tan 2 9i2- 

We shall see that a measurement of [p — p\ v - e to an accuracy better than ±3% 
is required in order to significantly improve our experimental knowledge of tan 2 612. 
The main reason for the required high precision in the measurement of [p — p] u - e is 
that existing experiments plus the luminosity constraint already determine the total 
p — p neutrino flux to ±2% (see the last column of table |J). 

We show in section [7.31 that a measurement of the p — p neutrino-electron scat- 
tering rate will provide an accurate measurement of the total luminosity of the Sun 
from neutrino measurements alone. 



In section [FTJ], we present for pep neutrinos the current predictions for the rate 
of v — e scattering and compare the power of measuring the pep scattering rate with 
the power of measuring the p — p scattering rate. 



7.1 Assuming the standard solar model p — p neutrino flux 

Table [?| shows what we can infer if we assume that the total p — p flux is known to 
an accuracy of ±1% from the standard solar model predictions. 

We show in the table the results of a global analysis if the expected v — e reduced 
scattering rate oip — p neutrinos, [p — p] u _ e , is measured to 5%, 3%, or 1% accuracy. 
The predicted rate for v — e scattering with a 50 keV electron energy threshold is 
(see table 2 of ref . @): 



[P-PU = 0.70 ±0.02(1°;°*). (7.2) 

We include in the global analysis that was used to construct table [7| all the presently 
known solar data, a simulation of three years of KamLAND measurements, and a 
measurement of the v — e scattering rate of 7 Be solar neutrinos, [ 7 Be] i/ c = 0.66(1 ± 
0.05). 



34 



Experiment 


s 


Am 2 


(io- 


~ 5 eV 2 ) 


tan 2 6> 


12 




fB 






fBe 




S + K 3 yr 


+ 7 Be 


7 q+0.2 
' -°-0.2 


K-0.6) 


0.41 


+0.04 , 
-0.04 \ 


+0.15\ 
-0.09J 


1.01 


+0.04 
-0.04 


f+0.13\ 
1-0.13 ) 


0.98 


+0.06 
-0.06 


'+0.19\ 
v-0.18j 


+ \P~P]v-e 


5% 


7.3^ 


-0.2 
-0.2 


1-0.6 J 


0.41 


+0.04 / 
-0.04 \ 


+0.14\ 
-0.09J 


1.01 


+0.04 
-0.04 


f+0.13\ 
\-Q.V2) 


0.98 


+0.06 
-0.06 


'+0.19\ 
.-0.18J 


+\p-p]u-e 


3% 


7.3^ 


-0.2 
-0.2 


K-0.6) 


0.42 


+0.03 , 
-0.03 I 


+0.12\ 
—0.11 J 


1.00 


+0.03 
-0.03 


V-O.llJ 


0.98 


+0.06 
-0.05 


'+0.19\ 
v-0.18j 


+ \P-P]u~e 


1% 


7.3^ 


-0.2 
-0.2 


( +0.7\ 
1-0.6J 


0.42 


+0.02 , 
-0.03 I 


+0.07\ 
-0.08J 


0.99 


+0.02 
-0.02 


(+0.07\ 
V-0.077 


0.99 


+0.05 
-0.06 


'+0.18\ 
v-0.18j 



Table 7: The effect of a measurement of the p — p v — e scattering rate: Assume 
standard solar model p — p neutrino flux. The table presents, in rows two-four, the 
best-fit and the la (3a) ranges computed by global fits that include a simulated measure- 
ment of the expected v — e scattering rate of p — p solar neutrinos. The p — p neutrino 
scattering rate is presumed to be measured with either 5%, 3% or 1% precision. In all cases, 
we have also included a simulated measurement with 5% precision of the v — e scattering 
rate of 7 Be neutrinos, a simulated result of three years of KamLAND data (K), and all 
available data from existing solar neutrino experiments (S) || 0, || [H], 16, [l?], |TJ, |33|] . The 
results from a global analysis that was made using just the solar, KamLAND, and 7 Be data 
are presented in the first row of the table. The p — p neutrino flux is constrained to have 
the same best-estimate value and uncertainty as in the BP00 solar model prediction |l4]] . 
The 8 B and 7 Be fluxes were treated as free parameters, but the luminosity constraint is not 
imposed. In constructing each column of the table, we have marginalized over all variables 
except the one whose range is shown. 



There are, according to table 0, only modest improvements in our knowledge of 
the neutrino oscillation parameters and the neutrino fluxes until the accuracy of the 
p — p measurement becomes better than ±3%. In other words, the measurement 
accuracy must be comparable to the quoted accuracy, ±1%, of the solar model 
calculation of the total p—p flux in order to make a major improvement. If [p — p\ v _ c 
is measured to an accuracy of 1%, then the 3a uncertainties in determining tan 2 612 
and /b will be reduced by a factor of two. The uncertainties for Am 2 and f-Qe will 
be affected only by relatively small amounts even if [p — p\ v _ c is measured to 1%. 

7.2 Free neutrino fluxes 

Suppose we treat the p—p, 7 Be, 8 B, and CNO solar neutrino fluxes as free parameters 
and impose the luminosity condition. Suppose also the p — p and 7 Be neutrino- 
electron scattering rates are measured. 

In this case, table [8] and figure |6| show that a ±1% measurement of the rate of 
v — e scattering rate by p — p solar neutrinos will improve by a factor of two our 
experimental knowledge of the allowed range of tan 2 612 , but will not significantly 
improve our knowledge of the other solar neutrino fluxes or of Am 2 . An accuracy of 
better than ±3% is necessary, according to table [§], to decrease the allowed range of 
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Figure 6: Allowed oscillation parameters: p—p + 7 Be solar neutrino experiments 
plus existing Solar neutrino experiments plus 3 years of KamLAND. The figure 
illustrates the effect of a future 3% (la) or 1% (la) measurement of the rate, ]p — p] u _ c , 
of p — p neutrino-electron scattering. The p — p, 7 Be, 8 B, and CNO neutrinos are treated 
as free parameters subject to the luminosity constraint [12]. Each contour is calculated 



making use of all the currently available solar neutrino data ||, 0, |8], [TJ, |TJ], [Tt], ^(], ||] plus 
the anticipated data from three years of operation of the KamLAND reactor experiment 
(cf. figure ||) plus a 5% measurement of the 7 Be solar neutrino flux (cf. figure |5|). 

tan 2 8 12 by more than 15%. 

A 1 % measurement of the p — p scattering rate will also improve the bound on 
the sterile neutrino fraction, assuming that there are no sterile neutrinos. Using 
the existing solar data and simulated three years of KamLAND data plus a 5% 
measurement of the 7 Be neutrino-electron scattering rate and a 1% measurement of 
the p — p neutrino-electron scattering rate, we find that 

sin 2 ?? < 0.08. (7.3) 



The result shown in equation [773] was obtained by allowing all of the solar neutrino 
fluxes to be free parameters and imposing the luminosity constraint. Thus a 1% p—p 
measurement can improve the limit on sin 2 rj by about 20% (cf. the discussion in 
section |4.4j) . A 3% measurement of the p — p neutrino-electron scattering rate will 
not significantly improve the limit on sin 2 rj. 

What can more precise solar neutrino experiments tell us about the mixing angle 
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Table 8: The effect of a measurement of the p — p v — e scattering rate: free 
fluxes, luminosity constraint. The table presents, in rows two-four, the best-fit and 
the la (3d) ranges computed by global fits that include a simulated measurement of the 
expected v — e scattering rate of p — p solar neutrinos. The p — p neutrino scattering rate 
is presumed to be measured with either 5%, 3% or 1 % precision. In all cases, we have 
also included a simulated measurement with 5 % precision of the v — e scattering rate of 
7 Be neutrinos, a simulated result of three years of KamLAND data (K), and all currently 
available data from existing solar neutrino experiments (S) || [7|, || [j~5|, [[(], 
The global analysis using just the solar, KamLAND, and 7 Be results is presented in the 
first row of the table. The p — p, 7 Be, 8 B, and CNO neutrino fluxes are all treated as free 
parameters subject to the luminosity constraint. In constructing each column of the table, 
we have marginalized over all variables except the one whose range is shown. The allowed 
range of Am 2 is the same for all cases shown: Am 2 = 7.3 ± 0.2 (log) eV 2 
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#13? A p — p solar neutrino experiment must be more accurate than 5% in order to 
improve the present limit that is derived from reactor, atmospheric, and solar neu- 
trino data (see figure |]and ref. [01). A p—p measurement accurate to 1% will, when 
analyzed together with a 5% 7 Be solar neutrino experiment, all presently available 
solar neutrino and reactor data, and three years of KamLAND data, be capable of 
constraining #13 to sin 2 9 13 < 0.036 at 3a. Figure |2| illustrates this modest improve- 
ment compared to the current limit of sin 2 ^ < 0.05. A 1% 7 Be solar neutrino 
experiment will not improve the bound shown in figure |^ without an accurate p — p 
experiment. 

We quote in all of our tables in this paper the allowed ranges after marginaliz- 
ing over every parameter except the one of special interest. We have, for example, 
marginalized over all other quantities except sin 2 77 in obtaining equation |7.3| . This 
marginalization procedure is not the common practice in the field of neutrino oscilla- 
tions. Because what we mean by errors is somewhat non-standard (but we think our 



procedure is valid), we have included in section [T4| and also in section |97J discussions 



of the effect of marginalizing over input parameters. 

7.3 How accurately will a p — p neutrino experiment determine the solar 
luminosity? 



We have discussed in section £05 and section |6.5| how neutrino experiments can be 
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used to measure the total energy generation rate of the Sun. We have carried out a 
global analysis of all of the existing solar data plus a simulated 5% measurement of 
the p — p neutrino-electron scattering rate as well as three years of simulated data of 
KamLAND and a 5% measurement of the Be neutrino-electron scattering rate. We 
find 



Lq (neutrino — inferred) 

Tq 



n qq+o.04 



(-ail) ( 5 % P-P experiment) . (7.4) 



Comparing equation |6.10| and equation [7.4| , we see that a 5% p — p measurement 
can improve the neutrino measurement of the solar luminosity by a factor of about 
three. 

The constraints become even a factor of two stronger if the p — p measurement 
is made to an accuracy of 1%. In this case, we find 



Lq (neutrino — inferred) n nn+0 02 / +0 07 



0-99±o:oa {-ool) (1% p-p experiment) . (7.5) 



A p — p measurement to an accuracy of 1% would make a truly fundamental con- 
tribution to our knowledge of stellar energy generation. The neutrino measurements 
could be used to bound to ±2% the contributions of all sources of energy other than 
the low-energy fusion reactions among light elements (the well-known p — p chain 
and CNO cycle). 

Unfortunately, a measurement of the p—p neutrino electron scattering rate to an 
accuracy of 1% will not improve the constraints on the fraction of the solar luminosity 
that is produced by CNO fusion reactions. That is to say, the limits summarized 
in equation |0] are not tightened, to the quoted accuracy, by a measurement of the 
p — p scattering rate to a precision of 1%. 

7.4 What about the pep line? 

In this subsection, we explore what can be learned by studying the pep solar neutrino 
flux. The interested reader can see in ref. [IS] an early discussion of the theoretical 



questions that can be addressed by combining the pep measurement with other solar 
neutrino measurements. Moreover, that has recently been renewed interest in the 
experimental possibilities for measuring the pep neutrino-electron scattering rate at 
a deep underground site (49]. 



We present in section [7.4. 1| the current best-estimate for the rate at which pep 



neutrinos are scattered by electrons. In section |7.4.2| , we compare the power of a pep 
measurement with the power of a p — p measurement. 

The ratio of the pep to the p — p neutrino flux is robustly determined by the 
standard solar model calculations. The ratio is determined more accurately than the 
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individual fluxes because the ratio only depends weakly on the solar model charac- 
teristics. We have therefore used the BPOO value for the ratio, 

~ g) = 425.9, (7.6) 

<p(pep) 



in the calculations reported in this subsection (and in section [7[T] and section |7.2|) 



However, we have verified that none of our quantitative conclusions are affected 



significantly if we change the ratio given in equation |7.6| by a relatively large amount, 
±10%. 

7.4.1 The current predictions 

What is the current best-estimate prediction for the rate of scattering of pep solar 
neutrinos by electrons? To answer this question, we have repeated the global analysis 



of the existing solar and KamLAND data as described in section |6J]. But in the 
present discussion we present the predictions for the pep, 1.4 MeV, neutrino line 
instead of the 7 Be neutrino line. 

If we treat the 8 B solar neutrino flux as a free parameter, but assume that 
all the other neutrino fluxes are as predicted by the BPOO standard solar model 



calculation [[L4}1 , then the current best-estimate prediction for the reduced pep reduced 



reaction rate is 

[pep] u _ e = 0.65±0.02(^). (7.7) 

If instead we treat all the neutrino fluxes as free parameters and impose the lumi- 
nosity constraint, then the predicted pep neutrino scattering scattering rate is 

[pep} u _ e = 0.66^(^ 4 7 ) [*H free] . (7.8) 
7.4.2 Which is better, p — p or pep? 

Assuming that the pep neutrino flux is measured instead of the p—p neutrino flux, we 
have repeated the global analyses of existing and future solar and KamLAND data 
that are described in section [7J] and section |7T2| . We assumed for these calculations 



the same accuracy for the measurements of the pep neutrino flux as was assumed in 



section TA and section T71 for the measurements of the p — p neutrino flux. 

Our global analyses show that a measurement of the v — e scattering rate by pep 
solar neutrinos would yield essentially equivalent information about neutrino oscil- 
lation parameters and solar neutrino fluxes as a measurement of the v — e scattering 
rate by p — p solar neutrinos. The estimated best-estimates and uncertainties in the 
parameters are almost identical for the analyses we have carried out for p — p and 
pep neutrinos. 



39 



8. Maximal Mixing? 



Is the currently allowed oscillation region for solar neutrinos consistent with maximal 
mixing, i.e., with tan 2 6* 12 = 1? And, if not, at what level is maximal mixing excluded? 
These are important questions for constructing particle physics models (see, e.g., 
[0, 50,1511153). 



We determine in section |871] at what level maximal mixing is excluded by using 
all the currently available solar plus KamLAND data. In section we investigate 



the extent to which the SNO salt-phase data |53], |54| will strengthen or weaken 
the conclusion that exact Maximal Mixing is disfavored by solar and reactor data. 
In order to illustrate the sensitivity (or insensitivity) of the results to individual 
experiments, we present the results of the global analyses with and without including 
the chlorine experiment. 

8.1 Solar plus KamLAND data today 

We see from figure ||] that the la allowed contour in the Am 2 -tan 2 # 12 plane satisfies 
tan 2 $i2 < 0.52. If we perform a global solution including all solar and KamLAND 
data and marginalize over Am 2 , then maximal mixing is excluded at 5.3a, i.e., 

tan 2 # 12 < 1 at 5.3a. (8.1) 



This result given in equation |S.1| is in good agreement with the value of 5Aa for 
the extent of the exclusion that was obtained [53J by the SNO collaboration using 
independent computer codes. 



The result stated in equation IO is slightly weakened if we exclude the chlorine 



experiment in the global solution. Then 

tan 2 #12 < 1 at 4.8a [without chlorine.] (8.2) 

If we had not marginalized over Am 2 , but instead had used the two-dimensional con- 
tours of Am 2 -tan 2 9\2 (as is sometimes done in the literature), then the results would 



have been weaker than the constraints given in equation |8.1| and in equation [O 
We would have obtained that, for 2 d.o.f. , maximal mixing was excluded at 4.9a 
(including chlorine) or 4.4a (excluding chlorine). 

8.2 How well did the simulations work? 



We described in section of hep-ph/0305159, v2. how we simulated, prior to their 
announcement, the results of the SNO measurements in the salt phase (see also the 
footnote in section [T6]of the present version of our paper). Using simulated instead 
of real SNO salt-phase data, we estimated that 

tan 2 #i2 < 1 at 5.5a [with simulated SNO salt — phase data]. (8-3) 
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Thus the simulations correctly predicted that the salt-phase results would strengthen 
the exclusion of maximal mixing by about an additional 2a relative to the exclusion 
of 3.5a that existed at the time the simulations were made. 



9. Discussion 



In section we answer the question of why do solar neutrino experiments below 1 
MeV. We use as the basis for our answer the calculations and discussion presented 
in the present paper. 

We summarize in section \9.2\ our view of what are the most important ideas and 
results that are described in section Q to section [|. In section |9l| we state our view 
of what all these results mean. 

9.1 Why do solar neutrino experiments below 1 MeV? 

There are three primary reasons for doing low energy solar neutrino experiments 4 . 

First, new phenomena may be revealed at low energies (< 1 MeV) that are 
not discernible at high energies (> 5 MeV). According to the currently accepted 
LMA oscillation solution, the basic oscillation mechanism switches somewhat below 
1 MeV from the MSW matter-dominated oscillations that prevail at high energies 
to the vacuum oscillations that dominate at low energies (More precisely, we mean 



by the phrase 'low energies' those energies for which (3 < cos26 , 12 , see equation \L3 
and figure [I]). We want to know if this transition from matter- induced to vacuum 
oscillations does indeed take place. If it does, is the ratio {(3) of the kinematic term 
in the Hamiltonian (i.e., Am 2 /2E) to the matter-induced term(v / 2G , F^ ) the only 
parameter that determines the physical processes that are observed in this energy 
range? Or, could there be entirely new physical phenomena that show up only at 
the low energies, very long baseline, and great sensitivity to matter effects provided 
by solar neutrino experiments (see, for example, refs. |56| . [571)? 

Second, new solar neutrino experiments will provide accurate measurements of 
the fluxes of the important p—p and 7 Be solar neutrino fluxes, which together amount 
to more than 98% of the total flux of solar neutrinos predicted by the standard solar 
model. These measurements will test the solar model predictions for the main energy- 
producing reactions, predictions that are more precise than for the higher-energy 
neutrinos. 



4 For earlier discussions of the reasons for doing low energy solar neutrino energy experiments, 
see ref. [j55| , e.g., and the many important talks in the LowNu conference series: 
http 
http 
http 

and http : //www. sns.ias.edu/ ~ jnb/Meetings/Lownu/ 



/ /cdfinfo.in2p3.fr/LowNu03/; 

//www.mpi — hd.mpg.de/nubis/wwwiownu2002/index.html; 
//www — sk.icrr.u — tokyo.ac.jp/lownu/; 
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Remarkably, the combination of a 7 Be solar neutrino experiment and a p — p 
solar neutrino experiment can determine experimentally which terminating reaction 
of the p — p chain, 3 He- 3 He or 3 He- 4 He, is faster in the solar interior and by how 
much. The ratio R of the rate of 3 He- 3 He reactions to the rate of 3 He- 4 He reactions 
averaged over the Sun can be expressed in terms of the p — p and 7 Be neutrino fluxes 
by the following simple relation 5 : 

< 3 He + 4 He> 20( 7 Be) 

< 3 He + 3 He > 0(pp) - 0( 7 Be) ' 1 ' ' 



The standard solar model predicts R = 0.174 [p8| . One of the reasons why it is 
so important to measure accurately the total p — p and 7 Be neutrino flux is in order 
to test this detailed prediction of standard solar models. The value of R reflects the 
competition between the two primary ways of terminating the p — p chain and hence 
is a critical probe of solar fusion. 

Using only the measurements of the solar neutrino fluxes, one can determine 
the current rate at which energy is being produced in the solar interior and can 
compare that energy generation rate with the observed photon luminosity emitted 
from the solar surface. This comparison will constitute a direct and accurate test of 
the fundamental idea that the Sun shines by nuclear reactions among light elements. 
Moreover, the neutrino flux measurements will test directly a general result of the 
standard solar model, namely, that the Sun is in a quasi-steady state in which the 
interior energy generation rate equals the surface radiation rate. 

Third, low-energy solar neutrino experiments, 7 Be plus p—p, will make possible a 
precise measurement of the vacuum mixing angle, #12, and an improved measurement 
of the sterile mixing parameter sin 2 r\ (as well as a slightly improved constraint on 
6*13) (see section |7.2| ). 

The low-energy measurements will provide, in addition, essential redundancy. 
In this paper, we have assumed the correctness of all solar neutrino and reactor 
experiments that have been performed so far or which will be performed in the 
future. But, the history of science teaches us that this is a dangerous assumption. 
Sometimes unrecognized systematic uncertainties can give misleading results. To 
be sure that our conclusions are robust, the same quantities must be measured in 
different ways. 

9.2 Principal results 

We present a list of our principal results, beginning with the results that pertain to 
individual solar neutrino fluxes (section |9.2.1| ), the neutrino oscillation parameters 
(section |9.2.2j ), and specific new experiments (section |9.2.3| ), and ending with tech- 
nical remarks concerning statistical significance and marginalization (section |9.2.4| ). 



3 More precisely, 0( 7 Be) should be replaced by the sum of the 7 Be and 8 B neutrino fluxes in the 



denominator of Eq. (9.1) 
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9.2.1 Solar neutrino fluxes 

The p — p solar neutrino flux. The existing solar and reactor experiments de- 
termine the flux of p — p solar neutrinos to an accuracy of ±2%, provided one 
imposes the luminosity constraint (see below and section |^). The measured 
value is 1.02 ± 0.02 times the flux predicted by the BP00 |L4] standard solar 
model (see the last two entries in the last column of table 

The T Be solar neutrino flux. The existing solar plus reactor experiments provide 
only loose constraints on the 7 Be solar neutrino flux. The 7 Be flux is restricted 
at la to the values 0.93^;|, corresponding to approximately a ±40% uncer- 
tainty (see the last row of table ^). The ±40% uncertainty results when one 
varies all the principal solar neutrino fluxes as free parameters subject to the 
luminosity constraint. 

The 8 B solar neutrino flux. The 8 B solar neutrino flux is robustly determined by 
the existing solar and KamLAND experiments. If one treats only the 8 B flux 
as a free parameter, then the experimentally-determined flux divided by the 
standard model (BP00) prediction is /b = 1.00 ± 0.04 (±0.13) (see the second 
row of table [J). The uncertainty is unchanged if one allows all of the neutrino 
fluxes to vary as free parameters, subject to the luminosity constraint (see the 
last row of table |]). 

The CNO solar neutrino fluxes. The existing solar neutrino experiments con- 
strain only very poorly the CNO fluxes. The 15 O flux is the most strongly 
constrained by the currently available data and even for this most favorable 
case the flux could be at la (3a) as large as 2.4 (6.4) times the standard model 
prediction (see equation |4.1| ). A 7 Be experiment accurate to ±5% could reduce 
the la (3a) uncertainty in the 15 neutrino flux by a factor of five (two) (see 
equation 



The power of the luminosity constraint. Table |4] shows that the luminosity con- 
straint reduces by a factor of ten the uncertainty in the experimentally-determined 
p — p solar neutrino flux, from ±0.22 (la) without the luminosity constraint to 
±0.02 with the constraint. One can understand from the approximate equa- 
tion ^TT] why the luminosity constraint provides such a powerful limitation on 
the p — p flux, but is not very restrictive for the 7 Be solar neutrino flux. If 
the real Sun is relatively close to the standard solar model about 85% 



of the energy generation in the model goes through the 3 He- 3 He termination 
(producing two p — p neutrinos) and only about 15% goes through the 3 He- 4 He 
termination (producing one p—p neutrino and one 7 Be neutrino). Therefore the 
condition for energy conservation that is expressed by the luminosity constraint 
limits the p — p neutrino flux more tightly than the 7 Be neutrino flux. 
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The luminosity constraint does not have a significant direct effect on the allowed 
region for the 8 B neutrino flux because the 8 B flux represents only a tiny 
fraction, 0.004%, of the energy generation in the standard solar model. The 
only effect of the luminosity constraint on the allowed 8 B range is indirectly 
through the allowed regions for tan 2 6u and Am 2 . 



Comparison of neutrino luminosity with photon luminosity. Neutrinos pro- 
duced by a particular nuclear reaction contribute an accurately known amount 
of energy to the Sun. Thus an energy-weighted linear combination of the indi- 
vidual solar neutrino fluxes represents the current rate at which energy is being 
generated in the solar interior [I2~|. From neutrino measurements alone, one 



can measure the solar energy generation rate and then compare this neutrino 
luminosity with the photon luminosity being radiated from the solar surface. 
This comparison would test the fundamental idea that nuclear fusion reactions 
are responsible for the energy radiated by the Sun. Moreover, this same com- 
parison would test a basic inference from the standard solar model, namely, 
that the Sun is in a quasi-steady state in which the energy currently radiated 
from the solar surface is currently balanced by the energy being produced by 
nuclear reactions in the solar interior. 

Using existing solar and reactor data, the neutrino luminosity can be deter- 
mined to about 20% (see equation |4.12|). A 7 Be solar neutrino experiment accu- 



rate to 5% could improve this determination to about 13% (see equation |6.10| ). 
The global combination of a 7 Be experiment, plus a p — p experiment, plus 
the existing solar data, and three years of KamLAND would make possible a 
precise determination of the solar neutrino luminosity. A p — p solar neutrino 
experiment accurate to 5% would make possible a measurement of the solar 



neutrino luminosity to 4% (see equation 7.4 ) and a 1% p — p experiment would 



determine the solar luminosity to 2% (see equation |7.5|). 



9.2.2 Neutrino parameters 

The allowed region of Am 2 and tan 2 12 . Figure ^ illustrates the currently al- 
lowed la region for Am 2 and tan 2 # 12 , as well as the much smaller allowed 
region that is expected to result after three years of operation of the Kam- 
LAND experiment. If only the 8 B neutrino flux is treated as a free parameter 
(all other neutrino fluxes and their uncertainties being taken from the standard 
solar model), then at lcr (3cr) (see the second row of table 0) Am 2 = 7.11q;4 
(± 2 ;e)and tan 2 8 12 = 0.42to;o4 (-oil)- If instead all the principal neutrino fluxes 
are treated as free variables but subject to the luminosity constraint, then (see 
the last row of table |) Am 2 = 7.3t^ (±11) and tan 2 12 = 0.4ll°;°l (±g;?|). 
The oscillation parameters are less well constrained if one requires that the 
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neutrino fluxes be determined by experiment rather than adopted from the 
standard solar model. 



A p — p solar neutrino experiment accurate to 1% can, when combined with 
existing solar and reactor data, a 5% 7 Be solar neutrino measurement, and three 
years of KamLAND data, lead to a modest improvement in the constraint on 
the mixing angle # 13 (see figure |2| and section [77 



The active-sterile admixture. The sterile fraction of the incident solar neutrino 
flux can be described in terms of the parameter sin 2 77, which is defined in 
section |3.2| and in ref. H. The la allowed range for the active-sterile admix- 



ture is sin 77 < 0.10 (see text before equation [4.9| ) if we allow all the neu- 
trino fluxes to vary as free parameters but impose the luminosity constraint. 
For the Super-Kamiokande and SNO experiments, this limit corresponds to 
/B.steriie = 0-0±nnn (see equation pjp. 

A 1% measurement of the p — p neutrino-electron scattering rate, in conjunc- 
tion with improved KamLAND measurements and a 5% 7 Be experiment, can 
significantly improve the existing limit on sin 2 r\ (see equation |7.3| ) provided 
that there are no sterile neutrinos that show up at low neutrino energies. The 
7 Be measurement determines the p — p flux significantly accurately that one is 
then sensitive to a small loss of flux to sterile neutrinos in a p — p experiment. 
However, even a 1% measurement of the 7 Be neutrino scattering rate will not, 
without an accurate p — p measurement, tighten significantly the bound on 
sin 2 rj (see discussion in section |6~2|) . 



The vacuum-matter transition. We define in equation [2.2| and equation [2.3| the 
quantity (3 [(3 = 2\^2G Fn e E u / Am 2 ], which governs the transition from vac- 
uum neutrino oscillations to matter dominated oscillations. For values of 
(3 < cos2# 12 , vacuum (kinematic) oscillations are dominant while for f3 > 1, 
matter (LMA) oscillations are dominant. 

Figure [l] illustrates the two limiting regimes, vacuum and matter oscillations, 
and the transition region between them. In the Sun, the vacuum-matter tran- 
sition occurs somewhere near 1 MeV. 

One of the principal goals of future solar neutrino experiments is to try to find 
direct evidence for the vacuum-matter transition. 



9.2.3 New experiments 

What additional KamLAND measurements will teach us. We have simulated 
the results of a full three years of KamLAND operations (see section |3.5| ), as- 
suming that the true result values of the oscillation parameters are equal to 
the current best-fit parameters. Performing a global analysis of the simulated 
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KamLAND 3 yr data and all the currently available solar neutrino data, we 
find that the additional KamLAND data will reduce the allowed 3a region for 
Am 2 by about a factor of two (compare the second and third rows of table |2] 
and table |3[). The allowed regions for tan 2 12 and /b will be reduced less 
dramatically, but still significant, by about a factor of one-fourth. 

A 7 Be solar neutrino experiment. The predicted rate (in terms of the standard 
model rate) for a 7 Be neutrino-electron experiment is very well determined if 
one assumes the standard solar model is correct: [ 7 Be] !/ _ e = 0.66 ±0.02 (la) 
(see equation |6.2|) . However, if one does not assume the correctness of the 



solar model, then the uncertainty in the prediction is about nine times worse: 
[ 7 Be] v _ c = 0.58ig;^ (see equation pD- 



We need an experiment to measure directly the flux of 7 Be solar neutrinos! 

How accurate does the 7 Be experiment have to be in order to provide important 
new information? This question is answered in table [l| A measurement of the 
u — e scattering rate accurate to ±10% or better will reduce by a factor of four 
the uncertainty in the measured 7 Be neutrino flux. Moreover, the 10% 7 Be flux 
measurement will reduce the uncertainty in the crucial p — p flux by a factor of 
about 2.5. A 7 Be measurement accurate to ±3% would provide another factor 
of two improvement in the accuracy of the 7 Be and p — p solar neutrino fluxes. 

All of these improvements are measured with respect to what we expect can be 
achieved with three years of operation of the KamLAND experiment (see the 
top row of table which is likely to be completed before a 7 Be solar neutrino 
experiment is completed. Comparable information can be obtained from a 
CC (neutrino absorption) experiment and from a neutrino-electron scattering 
experiment if both are performed to the same accuracy. 

Contrary to what some authors have stated, a 7 Be solar neutrino experiment 
is not expected to provide significantly more accurate values for the neutrino 
oscillation parameters than what we think will be available after three years of 
operation of KamLAND (see table |6|). 

A p — p solar neutrino experiment. According to the standard solar model, about 
91% of the total flux of the neutrinos from the Sun is in the form of the low 
energy (< 0.42 MeV) p — p neutrinos. We cannot be sure that we have an 
essentially correct description of the solar interior until this fundamental pre- 
diction is tested. Moreover, the p — p neutrinos are in the range where vacuum 
oscillations dominate over matter effects, so observing these low-energy neu- 
trinos is an opportunity to test in a crucial way also our understanding of the 
neutrino physics. 
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If we really know what we think we know, if the standard solar model is correct 
to the stated accuracy (±1% for the total p — p neutrino flux) and if there 
is no new physics that shows up below 0.4 MeV, then table |8] shows that a 
measurement of the p — p flux to an accuracy of better than ±3% is necessary 
in order to significantly improve our experimental knowledge of tan 2 8u. The 
main reason why such high accuracy is required is that the existing experiments, 
if they are all correct to their quoted accuracy, already determine the p—p solar 
neutrino flux to ±2%. (We assumed in constructing table [8] that three years 
of KamLAND reactor data will be available, as well as ±5% measurement 
of the 7 Be neutrino-electron scattering rate. The 7 Be measurement does not 
contribute significantly to the measurement accuracy for tan 2 6u.) 

As described above, an accurate measurement of the p — p solar neutrino flux 
will provide a direct test of the fundamental ideas underlying the standard solar 
model. The p—p measurement will make possible the determination of the total 
solar luminosity from just neutrino experiments alone. The neutrino luminosity 
can be compared with the photon luminosity to check whether nuclear fusion 
reactions among light elements is the only discernible source of solar energy 
and whether the Sun is in an approximate steady state in which the rate of 
interior energy generation equals the rate at which energy is radiated through 
the solar surface. 

The pep neutrinos (a 1.4 MeV neutrino line) can give essentially the same 
information as the p — p neutrinos. 

Maximal mixing; the SNO salt-phase data. The solar (without the SNO salt- 
phase data) plus reactor data disfavor neutrino oscillation solutions with max- 
imal mixing, i. e., tan 2 812 = 1, at a confidence level of 3.5<r. The addition of 
the most recent solar data, excludes maximal mixing at a confidence level of 
5.3cr (see equation ^.1|). 



What if one experiment is wrong? The implications of solar neutrino experi- 
ments for physics and astronomy are too important to be allowed to depend 
upon just one experiment, no matter how well that experiment appears to 
have been performed. We have therefore checked the sensitivity of some of our 
results to removing the chlorine experiment from the total data set. 

The principal reason for choosing the chlorine experiment to test the robustness 
of our conclusions is that the experiment supplies an important number. The 
total capture rate of solar neutrinos by 37 CI affects significantly some of our 
inferences. Although the chlorine experiment has been tested internally in 
different ways and has been carried out with exemplary care and skill JJ, 
the full experiment has not been directly checked by any other solar neutrino 
experiment. 
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If we include the chlorine experiment in the total data set, the reduced 7 Be neu- 
trino flux implied by all the data is /s e = 0.93^0 63 (see table §). If f ne chlorine 
experiment is removed from the data set then f Be = 1. 24ti%\ (see equation[4~2|). 
Thus removing the chlorine experiment degrades the existing 7 Be neutrino flux 
measurement from a ±47% measurement to a ±63% measurement. 

Removing the chlorine experiment from the data set results in a relatively 
slightly weaker conclusion regarding maximal mixing. Instead of maximal mix- 
ing being disfavored by 5.3a, as it is with the full data set, maximal mixing 
would only be disfavored by 4.8a if the chlorine experiment is omitted (cf. 
equation |8~J with equation p72| .) 



9.2.4 Statistical significance and marginalization 

Statistical significance. Whenever we present the allowed range of a particular 
variable [e.g., Am 2 , tan 2 12 , sin 2 9 13 , sin 2 ^, f B , f Be , f p _ p , or L (CNO)], we 
have marginalized over all parameters except the quantity singled out for special 



interest. This procedure is illustrated in section 3.3 , where the marginalization 



procedure is outlined specifically for the variable #13. The marginalization 



procedure has also been applied (see discussion in section |3.4| ) to other derived 
quantities like [ 7 Be]^_ e , or [p — p] v - e (or \pep] u _ e ) that are themselves functions 
of the fundamental input variables. 

The allowed ranges of [ 7 Be]„_ e or [p — p] v - e calculated by the marginalization 
procedure described above are different from, and smaller than (by up to about 
50%), the allowed ranges that were usually quoted in our own previous papers. 
There, we determined the uncertainties for quantities like [ 7 Be]^_ e or [p — p] v ~e 
by sampling the two-dimensional parameter space, Am 2 and tan 2 #12, at the de- 
sired confidence level, together with the other input parameters. The marginal- 
ization procedure described here yields a well defined range independent of the 
number of free parameters in the analysis. The marginalization procedure has 
been used to extract derived parameters as fB, s teriie @ or Am 2 sin 29i 2 (first ref. 
in 0) and it should be used in the future when calculating individual elements 
of the neutrino oscillation matrix like cos #13 x cos #12. 



Technical Details. We describe in section |^ the method we use to calculate the 
global x 2 , introduce the parameter that represents the active-sterile admixture, 
and define the reduced solar neutrino fluxes f B , fBe, and f p - p . We also summary 
the 96 data points (see table [I]) that we use in our global analyses. In addition, 
we outline how we create simulated data for three years of KamLAND operation 
and for the salt-phase measurements of SNO. 
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9.3 What does it all mean? 

We are agnostics on the question of whether or not we know the physics and the as- 
tronomy well enough to extrapolate accurately to a different neutrino energy domain 
and to different solar neutrino sources. New physics or astronomical processes may 
be relevant at neutrino energies below 1 MeV that are not important at higher ener- 
gies. We have made the extrapolation in this paper in order to have a quantitative 
basis for discussing what to expect and what we might learn. 
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